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Abstract 



The energy spectrum of the baryons is determined by treating each of them as a three-body system with the 
interacting forces coming from a set of two-body potentials that depend on both the distance between the 
quarks and the spin and orbital angular momentum coupling terms. Constraint dynamics is first reviewed 
for a relativistic two-body system in order to assemble the necessary two body framework for the three-body 
problem and then we review the different types of covariant two-body interactions involved in constraint 
dynamics, including vector and scalar, and solve the problem of energy eigenstates using constraint dynamics. 
The Two Body Dirac equations of constraint dynamics derived by Crater and Van Alstine, matched with the 
quasipotential formalism of Todorov as the underlying two-body formalism are used, as well as the three- 
body constraint formalism of Sazdjian to integrate the three two-body equations into a single relativistically 
covariant three body equation for the bound state energies. The results arc analyzed and compared to 
experiment using a best fit method and several different algorithms, including a gradient approach, and 
Monte Carlo method. 
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1 Introduction 



Recent quark model calculations done by Crater et al. [H 121 13] using covariant two-body Dirac equations 
in a relativistic constraint dynamics formalism have given a good description of the meson masses for both 
light and heavy quarks, using world scalar and vector potentials that depend on merely one or two param- 
eters. The good quality of the fit has been attributed to the exact two-body kinematics merged with a 
QCD interaction potential based on vector and scalar potentials that uses a minimal number of variable 
parameters. These non-perturbative, i.e. numerical results hold up well when compared to other methods 
for meson spectroscopy. The vector potentials, in turn, have a structure originally derived from the classical 
electrodynamics of Wheeler and Feynman.[4 This structure can also be obtained from quantum electro- 
dynamics by using a covariant three dimensional truncation of the Bethe-Salpeter equation based on the 
Todorov quasipotential, which is then compared to the Two-Body Dirac equations. [3] The comparison is 
done in order to identify the appropriate invariant potential functions that will be used in the potential 
model. 

A usual common ground for a two-body bound state problem (almost regardless of which approach is 
taken) is the Bethe-Salpeter equation, which is generally not used in its full (in other words, four dimensional) 
form due to difficulties with the relative time coordinate. The Two-Body Dirac equations of constraint 
dynamics give a covariant, three-dimensional version of the Bethe-Salpeter equation that solves the problem 
of relative time and relative coordinates. Originally, the Two-Body Dirac equations of constraint dynamics 
arose from a supersymmetric treatment of two pseudoclassical constraints (with Grassmann variables in place 
of gamma matrices) for scalar interactions, which were then quantized. Later work led to the development of 
the form of the Two-Body Dirac equations that were not limited to scalar interactions. Covariant constraint 
dynamics, at its core, is a method which allows connections between non-independent operators to be built 
in a system that has more variables than it has observables, such as occurs with the mass shell constraint 
+ ~ 0. It is inherently four-dimensional, but only three of the momentum variables are independent. 

a 

In this work, we review the Hamiltonian constraint dynamics formalism for the two-body system and 
apply it to the three-body quark problem for baryon spectroscopy. The derivations of the Two-Body Dirac 
equations are reviewed along with their potentials for a relativistic, spin-dependent formalism. In taking 
the two-body equations to a three-body system we still regard the system as the naive quark model in 
that each interaction is between each pair of quarks and there is no over-arching three-body interaction 
to be considered, but the system now has three sets of interactions instead of just one. Thus, all of our 
interactions are still two-body interactions, but for three sets of quarks. This works out relatively simply for 
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some of the interacting potential but has required extensive reworking for others, such as the spin-spin and 
spin-orbit interactions. Not surprisingly then, a large portion of this work is devoted to the development of 
the two-body theory as this provides both the tools for working the three-body case along with the possible 
drawbacks to consider as we the formalism is extended to larger systems (e.g. compatibility of the constraint 
equations). The end result of this, as will be seen later, is that the theoretical considerations that must 
be applied in order to go to the three-body system are relatively simple, but the actual mechanics of doing 
so are not and multiple methods that are concise and rigorous (such as separation of variables and relative 
coordinates) in the two-body case cannot be applied to the three-body one, at least not without a substantial 
adaptation of the two-body formalism. 

The Hamiltonian Constraint Dynamics formalism 'W allows for a relativistic method of accounting for 
two-body effects. In addition. Constraint Dynamics as developed by Crater and Van Alstine 4 provides 
not only the usual spin interaction dependence seen in the Dirac equation but also additional terms needed 
to make the approach mathematically consistent. It is correct and useful in both a classical and a quantum 
mechanical formalism and as such we show how it is derived and applied for a two-body system. Furthermore, 
we show how the world scalar and vector interactions come about in this formalism and how they are used. 
All of the interacting potentials that are used in this work (and the two-body system as well) are dependent 
on these vector and scalar potentials. 

The variational principle is used with a Gaussian basis wavefunction of total JM to solve our eigenvalue 
equation. We show how the Gaussian basis can be expanded in a variational theory in order to achieve an 
eigenvalue closer to the true ground state. The matrix must be truncated after a reasonable limit is reached, 
but in theory an infinitely large variational matrix would give the exact ground state. Since the potentials 
are dependent on the energy eigenvalue, the standard approach to using the variational principle does not 
readily apply. A recursion algorithm is therefore used so that the effective Hamiltonian H has an embedded 
E dependence, but this E will change as we approach convergence. The program then is designed to 
iteratively solve these equations until a desired level of convergence is reached. As in the two-body case, the 
solution to our Hamiltonian can be done exactly for the kinematics but including the interacting potentials 
requires a numerical treatment. Also as in the two-body case, these interacting potentials depend only on 
one or two parameters. The numerical fitting routine uses a chi-squared minimization Monte Carlo routine 
combined with a simplified gradient approach to acquire a best fit for the spectrum of known baryons. We 
have compared our numerical results to both experimental data and to other theories, most notably the 
approach of Capstick and Isgur([B]), while also comparing the quark masses and potential parameters we 
obtained in our fit with those found in two-body results by Crater et al. [I] [2 [3] . Ideally, since this model 
uses the same basic model for the three-body system as has been used for the two-body system, one would 
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expect the parameters to be similar. The actual baryon masses are, of course, not compared to the meson 
results, but the masses of the quarks and values for other shared parameters are. 

1.1 Quarks and the Quark Model for the Baryons 

This section contains a brief review of relevant definitions of quarks and baryons as well as references 
to recent work on the baryon three-body problem. Quarks were introduced in the beginning of 1964 by 
Gcll-Mann and Zweig[7] as an idea evolving from considerations of unitary symmetry. Quarks are considered 
to be the fundamental building blocks of hadronic matter and each hadron is a bound state of quarks and/or 
antiquarks. Quarks are considered to be fermions in that they have odd half-integer spins. 

The baryons are three-quark bound states having an odd half-integer total angular momentum j. They 
experience the strong nuclear force, and are described by Fermi-Dirac statistics. The lowest-lying baryon 
states are three-quark states with orbital angular momentum / = 0, of which there are 8 states with total 
j — \ (the octet). This work begins with a focus on the octet, but quickly proceeds to a much larger number 
of baryons, such as the decuplet, which consists of Z = three-quark states with total j of | [7 , as well as 
most of the other baryon radial and orbital excitations (see Appendix B for the complete listing). 

Quantum field theory, from which the Bethe-Salpeter equation is directly derived, originated in the 1920s 
from the problem of combining a quantum mechanical theory with an electromagnetic field [8 . The theory 
was constructed by expressing the field's internal degrees of freedom as an infinite set of harmonic oscillators 
and applying a canonical quantization to them. The theory, of course, needed to combine relativity and 
quantum mechanics for the electrons as well, which is the reason for the importance of the Dirac equation. 

Quantum chromodynamics (QCD), which plays an important role in this work, is the theory of the strong 
interaction. The strong interaction is the force responsible for describing the interactions of quarks and 
gluons inside of a hadron and QCD itself possesses the properties of both asymptotic freedom and linear 
confinement. [8] Asymptotic freedom, which in our model comes from the vector potential term, allows for 
quarks and gluons to interact very weakly in high energy interactions where the force is approximately 
independent of r. However, quarks also exhibit linear confinement and it would take an infinite amount of 
energy to separate two quarks, which is modeled in our scalar potential term. 

As far as work done with 3-body systems similar to this work, Capstick and Isgur have studied the 
three-quark system in a relativized quark potential model with chromodynamics. Toward that end, they 
have developed a model that supports the phenomenology of nonrelativistic calculations. It excludes the 
spin-orbit interactions in baryons and leads to spectra and internal compositions similar to those of the 
nonrelativistic model ^ . In a similar fashion to our work, they use a set of two-body potentials that have 
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been previously derived and tested for the meson spectrum, as we do here with the meson model created 
by Crater et al.([TJ [31 [3]). It still does follow similar qualitative (for spectra and internal compositions) 
properties as the nonrelativistic quark model, since the goal of that work was to overcome the shortfalls 
of the nonrelativistic approach, such as the nonrelativistic treatment of both the quark motion and quark 
dynamics. Isgur and Karl have also studied the ground state baryon octet in a quark model with flavor- 
independent confinement and color hyperfine interactions in an even earlier study using a similar model to 
the later model of Capstick and Isgur, with good results for the octet [9J . Neither of these methods, however, 
use the relativistic constraint dynamics formalism used here, which allows us to create a completely covariant 
formalism. Capstick and Isgur also use basis functions that are always antisymmetric under exchange of 
quarks 1 and 2, which is a simplifying approximation we do not use, though it is worth noting it is useful 
for the vast majority of the spectrum (most have 2 quarks that are the same flavor). 

Loring et. al. have written a series of three papers treating light baryon resonances within a relativistically 
covariant quark model based on the Bethe-Salpeter equation and demonstrate how to solve the Bethe- 
Salpeter by reduction to the Salpeter equation [10] . Again, this is a significantly different formalism from 
our approach, but as we stated above, the Bethe-Salpeter equation tends to be a common point in most 
approaches to this problem. 

By contrast to those approaches, Sazdjian uses the manifestly covariant formalism with constraints for 
the construction of relativistic wave equations which describe the dynamics of N interacting spin or spin 
i particles with N > 2. The system is governed by a single dynamical wave equation that determines the 
eigenvalue of the total mass squared of the system. For N = 2 his equations are essentially the same as 
presented below and his notation has been adapted accordingly. The interactions in the general N case are 
introduced through two-body potentials, but many body potentials can also be incorporated Sazdjian's 
N body eigenvalue equations provide a useful starting point for work on the manifestly covariant relativistic 
three-body problem. 

In this work, some of the simplifying substitutions derived by Isgur have been used [9] and the interaction 
potentials have been treated as separate two-body interactions. The constraint dynamics formalism that 
has been successful for the two-body problem applications to meson spectroscopy is reviewed and used, |12j 
with a focus on how to extend this to the relativistic three-body formalism for the baryons. In addition 
to not assuming that two quarks are the same mass as Isgur does, we also do not consider full three-body 
interactions, choosing instead to split our problem into just three two-body interactions. 
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1.2 History of Dirac's Equation and its Two-Body Counterpart 

As with any good discussion of things quantum, we begin with the non-relativistic Schrodinger equation 
for a single particle 

^^{r, t) + V^WM'Cr, t) = ih^/^ir, t). (1) 

This equation, developed in 1926, forms the foundation for modern quantum mechanics, where p = —ih^. 
Natural units will be used, setting h — c — 1. This equation is not relativistically covariant (meaning that 
it does not treat r and t equally), but there are very well-established methods to solve this equation. Thus, 
if whichever wave equation is used can be modified to fit this form, then the solution is essentially known. 
The ability to put an equation into a Schrodinger-like form will prove to be invaluable in solving far more 
complicated systems. 

Since the Schrodinger equation is non-relativistic, the next obvious step comes in relativizing it. By 
taking the standard relativistic equation for total energy, E'^ — p^c^ + m^c^ = P^ + m^(in natural units), 
and substituting the quantum analogs of p and E, one arrives at the Klein-Gordon equation for particles 
without spin [14] 

(p2+m2)* = 0, (2) 

^ where = — (V^ — ^)(the d'Alembert operator). This equation looks promising in that it treats x 
and t symmetrically, however, it has difficulties with negative energy states and the interaction form did not 
include a spin term, which is a necessity for this work. 

After seeing the issues with the Klein-Gordon equation, Dirac then set himself upon fixing the problems 
it had, namely the negative probability densities due to the lack of first-order derivatives in time, not to 
mention negative energies. This lead to the development of his equation in 1928, which is 

(7>M + H* = 0, (3) 

where 7'^ are the Dirac matrices. This equation is fully relativistic, consistent with the principles of 
quantum mechanics (including positive probability densities), and designed for spin-^[15^ particles. It does 
not, however, remove the issue of negative energy solutions, but it does provide an explanation for them in 
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quantum field theory, in which the solution for a negative energy electron is interpreted as a solution for 
a positive energy positron. Thus, it is a first-order equation in time (like the Schrodinger equation) and 
accounts for the existence of antiparticles. In the presence of world scalar (S) and four- vector interactions 
{A), these equations become 

{{p - A)2 + (m + Sf) ^ = (7^ (p - A)^ + (m + 5)) ^Z- = . (4) 

" ' V ' 

Kloin-Gordon pi^.^^ 

The next logical step is then to create a two-body equation for spin-i particles since there now existed 
a one-body equation. Breit was the first to develop a 2-body Dirac equation in 1929, formed by summing 
two free-particle Dirac Hamiltonians with a spin dependent interaction obtained by substituting Dirac a 
matrices for velocities in the semi-relativistic electrodynamic interactions of Darwin fTBl 1171 I18j . given by 

a 1 

= {oLi ■ Pi -)- P-^mi a2 ■ P2 + /32"^2 [1 - 7:(ai • 0:2 + ■ f "2 ■ ?)]}* • (5) 

r 2 

While this equation was useful for the perturbative calculation of the electromagnetic bound-state spectrum 
of multi-electron atoms, which is what Breit designed it for, it was not a well-defined wave equation, meaning 
(for terms beyond Coulomb) it could only be used perturbatively, and it was also not relativistically covariant. 
In an attempt to reproduce results from the Breit equation in weak-coupling perturbation theory, Bethe 
and Salpeter derived a two-body wave eg nation [19] from quantum field theory. While their equation is 
in agreement with quantum electrodynamics, its use as a relativistic bound-state equation has difficulties. 
Nakanishi|20| discovered that the Bethe-Salpeter equation would give negative-norm solutions as a result of 
relative time, a relativistic degree of freedom. 

Todorov created an effective potential method for the two-body system [211 [12] by generalizing the Ein- 
stein condition i?^ + to an effective particle of relative motion. He attempted to solve the problem 
that the Bethe-Salpeter equation had with relative energy and relative time by creating an equation for a 
particle of relative motion relativistically, but three-dimensional. He introduced the relative momentum, 
an energy, and a mass of the particle of relative motion that satisfied the relation of the Einstein condition 
between mass, energy, and momentum. 

In order to address the problem of relative time and energy, Todorov, Kalb, Van Alstine, and others [25] 
derived covariant, canonical mechanics for the two-body system from Dirac's Hamiltonian constraint dynam- 
ics. In such a system the dynamics are given by a set of constraints that act both to constrain the motion 
in phase-space and to generate the Hamiltonian. The development of Dirac's constrained Hamiltonian for- 
malism for the two-body problem to include spin with scalar, vector, tensor, pseudoscalar, and pseudovector 
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interactions for spin - i particles has been established by Crater and Van Alstine[231 HSl 111 HSl ISZl [2H] and 
used extensively in this work. 



2 Hamiltonian Constraint Dynamics 

This section begins our review of the Hamiltonian constraint dynamics formalism. We begin with 
the two-body system in order to fully develop the set of tools that are used in the three-body system, 
but it will become apparent that the two-body methods require extensive reworking in order to use the 
same basic formalism for the three-body problem. This formalism is useful because it provides a method 
of accounting for two-body relativistic effects that incorporate the automatic spin interaction dependence 
that the Dirac methodology provides and also provides a covariant, three dimensional truncation of the 
Bethe-Salpeter equation. We review constraint dynamics first for a classical system and then for a quantum 
mechanical system. We derive the constraints that arise from mass shell conditions and provide a review of 
the interacting potentials. We first derive the entire system of constraints and interacting potentials for a 
classical system and then we show how it is quantized. In this section we describe how the constraints arise 
in the context of relativistic classical mechanics. The Lagrangian for a free particle 



L — — m-\/— i^, (6) 



has canonical four-momentum 



dL mx , , 

ox V-a;^ 



and Legendre Hamiltonian 



Hl = P^i"" -L= -j= + = 0. (8) 



From = "^^i = —rn?, one arrives at the mass shell constraint with {p =(e, p)) 



'Ha=P^ + m^ = -e^ + « 0, (9) 

with the weak equality sign w indicating that the equality can only be imposed once equations of motion have 
been evaluated. When one has constraints, the Hamiltonian, called the Dirac Hamiltonian, is 'Hd = 'Hl + 
Ax constraints. Note that when the Legendre Hamiltonian vanishes, the constraint itself serves as the 
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Hamiltonian, % =p^ + m^. In this context we can introduce interactions by using a generalized mass shell 
constraint 

p2 + m2 + $w0, (10) 

where $ is some interacting potential to be determined later. These equations provide the necessary 
definitions for discussing the two-body Hamiltonians which form the basis for most of our work. 



2.1 Two-Body Hamiltonians for Relativistic Constraint Dynamics 

In order to understand the relativistic three-body problem one must first be very familiar with the 
relativistic two-body one, since the three-body formalism we use is based on two-body interactions (our 
three-body problem is essentially written as three two-body problems) . This section therefore describes the 
methods used to set up two-body Hamiltonians. In the two-body case, in analogy with what occurs in the 
1-body case, we start with a 2-body Lagrangian that would lead not to one but to two mass shell constraints: 



H20 = pl + mj^O, (11) 
and as in the one-body case interactions would be introduced by mass shell constraints of the form 

■Hi = p? + m2 + $i(x,P) wO, 

■H2 = p^ + m2 + $2(a;,P)«0, (12) 
where x is the relative coordinate and P is the total momentum, defined as 

X = xi — X2, (13) 

P = P1+P2- 



The invariants $i(a;,P) are called quasipotentials due to having center of momentum (cm.) energy 
dependent potentials that describe deviations from the free mass shell constraint. Here again, the weak 
equality sign means that these constraints are only to be applied after the Poisson brackets (see e.g. Eq.pS]) 
below) have been evaluated. The Hamiltonian in terms of the constraints only is [29] 
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H = XiHi + X2H2, (14) 

where Aj are unknown Lagrange multipliers. The constraints are constant in time, meaning that {Hi, %} = 0, 
provided that 

{Ur,U2}^Q, (15) 

which is known as the compatibility condition. If we explicitly work out these Poisson brackets, we arrive 
at the equation 

2pi • {pi, $2} + 2pi • {$i,P2} + {$1, $2} « 0. (16) 

One assumes that the invariant functions are dependent on two coordinates and the total energy of the 
system, as given by 

^i = *i(^,^,«^), (17) 

in which w is the total cm. energy, 



P 



P 

w 



= -1, 



X± = X — X\\^ 

P-x± = 0, 



{a;^,p''} = {jjf''' - P^'P^') . (18) 

Here, a:|| and x± are the components parallel and perpendicular to the total momentum, respectively. Thus, 
the compatibility condition becomes 

9$2 d^2 d^i 

-^Pi-x±-^- 4pi • a^ii^ - 4p2 .x^-^-Ap2- 2^11^ + {$1, $2} « 0, (19) 
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for which the sohition is (in the simplest case) 



2 

$1 = $2 = 



which gives 



{^1,%} = -4P • = 0. 



Since the non-rclativistic relative momentum is defined by 



m2Pi - miP2 mim2 , , 

P = : = ; (vi - V2) = /XV , 

nil + 1^2 mi + 7712 



it is natural to take as its relativistic counterpart 



— {£2Pi - £iP2) , 
w 

where ei and £2 are the center of momentum energies of particles 1 and 2, respectively, defined 



El + £2 = W. 



Then Eq.^ and Eq.^ imply 



Pi = [-p^eiF+p and p2 = P = S2F-p, 

w w 



which leads to 



2 

'Hi-'H2^pl-pl + ml-ml^ {el - ^1) + " (^1 + £2) -P • P + - w 
Since the £i are the cm. energies, this implies 

-pi-P -~Pi~{pi-P2) -P^ -P1+P2 w'^+mj-ml 



£1 = 
£2 = 



WW 2w 2w 

-P2-P -pi - {P2-Pi) -P^ -pI+pI w'^+ml-ml 



WW 2w 2w 

which in turn gives 



£1 - £2 



9 2 
nil ~ "^2 
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If instead one defines 



2 2 
TOT — 

81-62^^ ^, 29 

w 



tfien 

e\-el = ml- ml (30) 
and (|26l) is a derived constrainiy with the simple form 

P-pwO. (31) 

For the remaining constraint one may use any combination 

H=AiHi+A2H2, (32) 

of the constraints. The same constraint is obtained provided that Ai + A2 — 1. The reason is that the 
constraint "Hi — ^2 = 2F • p « 0. Thus, 

Ui = pl+ni\ + <^ = 2eip-P + p^ -e\+m\ + <^ (33) 
« p2 - 6^ + $ 

« ^2 = -2e2p-P+/-e2+m2 + $ 

where 

=el-ml=£l~ml^ X{w^ ,ml,Tnl)/Aw'^ , (34) 
with the triangle function A is defined as 

\r 2 2 2\ 4, 4, 4 o22 o22 o22 rot:\ 

A(w , nil, 1712) ~ w + nil + ™2 ~ ^""^ "^1 ~ ^""^ "^2 " 2TO]^TO2. (35) 

Thus we have 

Hi^H2~H ^p^ - 6^ + (36) 
as (weakly) equivalent constraint forms. Without interactions, the solution to Eg. ([55)1 becomes 



^The condition Eq. ||20|I is known as the relativistic version of Newton's third law. As stated there it is a strong equahty. 
It can be relaxed to the weak form $1 — $2 = P ■ Px{^±) ~ ^^'^ still lead to compatible constraints. 
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w = \Jp'^ +m\ + ^Jp^ +m\, (37) 

indicating exact relativistic two-body kinematics in the center of momentum system (in the cm., frame, 
p = (0, p)). Thus, we now have a classical two-body Hamiltonian with exact relativistic two-body kinematics 
and interacting potentials, though the potentials are still in a general form. Therefore, we must now 
determine the specific form of the interacting potentials $, first for world scalar and vector interactions. 



2.2 World Scalar Interactions 

This section discusses the world scalar interaction, beginning with a single, spinless particle[25 . We 
first derive the scalar interaction for a one-body Hamiltonian and then use that result to derive the two-body 
analog of said Hamiltonian. In addition, we use knowledge of the Klein-Gordon equation (and thus what 
the equations must reduce down to given certain conditions) to derive the exact form in which the scalar 
potential appears. The substitution m — > m -f S, defined as Af, is used for a particle in an external scalar 
field. The Lagrangian, the momentum, and Legendre Hamiltonian are then 



L = -{m + S)y/^, (38) 



Mi 



Mx 



2 



Ul = + AfV^ = 0. (40) 

V — 

The momentum equation then directly leads to the constraint for the one-body case 



p2 + ^ (m + S')2 « 0. (41) 

For the two-body system, one begins with the free constraints of Eq.(ITT|). The fact that these constraints 
need to reduce down to a Klein-Gordon form in the case of one mass being much greater than the other 
mass implies that scalar interactions should be introduced by rrii — )■ Mi, 
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where 



H20 ■H2^pl + M^^pl + ml + <^2{.x^,P)^Q. (42) 

M, = TOi + S'i. (43) 



From the compatibihty condition Eq. (|T5|) and Eq.(|20 



the solution of which is 



$1 
$1 
$2 



$2 = 



2miSi + Si , 



M| - m2 = 2m2S'2 



(44) 



resulting in 



/r2 



(45) 



Ml = mi cosh £ + 7712 sinh £ , (46) 

M2 — 7772 cosh C + 777i sluh £ , 

where C, S'i,and S2 are functions of xj_ in terms of a single underlying invariant function governing the scalar 
interaction which will be called S{xj_). 

In order to determine C, Si, S2, and thus $s one constructs a mass shell condition for a single particle 
of relative motion analogous to 

+ {m + Sf mO. (47) 
This single-particle Klein-Gordon equation has the form 

p2 + $ ^ £2 _ ^2^ (48) 
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where 



$ = 2mS + S^. 



(49) 



To continue further, there need to be analogous equations for the two-body case. For a two-body system 
one has (from Eg. ([55)) ) 

p2 + $ « b^, (50) 

To parallel what is found in the one-body Klein-Gordon case, one uses the energy and mass for the effective 
particle of relative motion introduced by Todorov. These are defined by 



which satisfy the Einstein condition 



Thus 



W — mj — 1712 

mim2 
w 



b'^el-ml. (52) 



which is directly analogous to the one-body case Eq.(j 

Let us define the four vector momentum analogous to the single particle p'^ 



In the cm. system this is 



(51) 



+ ^ ei-mi, (53) 



pA" = _^ fliipM . (54) 
w 



P''«(0,P) + — (u;,0) = (£„,p) . (55) 
w 
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Since 



2 2 



the Hamiltonian Ea. (P5)) 



e^-m„, (56) 



n=p^ -b^ + <Ps ~0, (57) 



then assumes the Klein- Gordon form 

■H = p2 + m^ + $s. (58) 

By looking at the single particle Klein-Gordon equations Ea. (P5| . Ea. p7| . and Ea.(P^ one can see that the 
proper choice of $s to bring it to the form required is 

$s = 2m„5 + 52, (59) 

which finally gives us a Hamiltonian of 

n = p^ + (m^ + Sf =p^ -el + [m^ + Sf « 0. (60) 

Note that this constraint can be obtained from its free form "Ho = + m^"^ by replacing with m^, + S. 
This structure for the interaction is supported by two independent field-theoretic arguments, one classical 
and one quantum-mechanical|30l [31 I32j . In addition, the forms given in Eq. (j5ip are confirmed in both 
classical and quantum field theor jifj. [48] This formalism is still classical, but we will derive the quantum 
analog of it and the other potentials. 



2.3 World Vector Interactions 

This section describes the methods of treating the world vector {A^) interaction, which will include 
time-like and electromagnetic-like vector interactions as well. In a similar manner to our treatment of the 
scalar potential, we derive a modified mass-shell constraint with a vanishing Legendre Hamiltonian. Unlike 

^This 111! structure and these kinematical forms for the effective mass and energy arise from eikonal approximate sums of 
ladder and cross ladder diagrams arising from a scalar field theory. I28| 
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the scalar potential, there are multiple types of vector interactions and we will review the time-like and 
electromagnetic-like vector interactions in addition to deriving the basic form of the potential. Lastly, we 
will combine the electromagnetic-like vector interaction with the scalar interaction. 
With the choice [1 [55] 



L — —nwj—'jp' -|- A-x , (61) 
one defines the interaction for a single particle with a vector field. The conjugate momentum is 

dL mx , 

ox y/-X^ 

which is rewritten as 



7r = p-A=-^==. (63) 



This equation leads us to the mass shell constraint of 



n^{p- Af ^0 , (64) 



which from Eq. (|63|) becomes 



•H = (tt^ + m'^) w 0, (65) 
and now, as before, the Legendre Hamiltonian is zero. 

. 9 

px — L = _ — Ax — mv — -|- Ax = . (66) 

Similarly to what was done for the two-body scalar case, one writes in four vector notation the classical 
minimal substitution 

pf^pf-^r^'Tf, (67) 

which causes the free particle constraints to become the interacting ones below 
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Ui = 

■H2 = 



(68) 



Assuming translational invariance so that the potentials depend on coordinates only through the relative 
coordinate x, these have the form 



■Hi = p2_^m2 + $i(a;,pi,p2)«0 , (69) 
H2 = p^ + m2 + $2(a;,Pi,P2)«0 . (70) 

Now that the basic form has been derived, we continue to the specific time-like and electromagnetic-like 
vector interactions. 

2.3.1 Time-like Vector Interactions 

For time- like vector interactions one sets[ll [53] 



A = i?2P^-p^ (71) 
which, using the above definition (Eg. lp^ ) for the Hamiltonian with a vector potential, leads to 

Hi = TTl + ml=pl + ml + ^iK 0. 

•H2 = TTl + ml=pl + ml + ^2^Q- (72) 

where 
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$1 = Trl-pl = -Ef + sl + 2Eip-P = -2pi-Ai+Al^2eiVi-V'}, 
$2 = nl~pl^-E^+el-2E2p-P^~2p2-A2+Al^2e2V2-Vi. 



(73) 



Note that these are the vector equivalents to Ea. (|42| . Ea.(|43)). and Eq. (|44)) . respectively. 
The compatibihty condition, here in weak form (see above footnote), is 

$1 « -El + £2 « $2 ~ -Ei + el = $v(2:±), (74) 

where '^v{x±) is defined by 

^v{x±) = -El + el^ -El + el, (75) 



Et-E'2^ei-ei (76) 



The solution for Eq.(l76l) is 



El = £i cosh J + 62 sinh J , (77) 
E2 — £2 cosh J -\- El sinh J , 

where , Vi, and V2 are functions of xj^ in terms of a single underlying invariant function governing the 
scalar interaction which will be called V(x_l). 

Now we want to find the analogy of Eq. (|5D)) for time- like vector interactions. The one-body Klein-Gordon 
equation for this is 

-{E- Vf + « 0, (78) 

which can also be written as 
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with 



$ = 2EV - 



(80) 



The two-body constraint of Eg. ([55]) 



ni^n2~n=p^ -b^ + $v(-t±), 



(81) 



assumes a Klein- Gordon form 



n = p'^ + ml + <i>v = p'^ - el + ml + <i>v ^ 0. 



(82) 



By looking at Eq.^, Eq.^, and Eq.dHOl) one can see that 



n 



p2 + ml + 2e„V -V^=p^+ml 



W 



V)^ «0. 



(83) 



Note that this constraint can be obtained from its free form "Ho = + m^^ by replacing p with p — VP or 
Ew with £.uj — V. 

2.3.2 Electromagnetic-like Interactions 

For spinless particles, there are four independent vectors : xi , X2 , pi , and p2 ■ This means the A'^ can be 
written as [26] 



Four vector terms that would produce unobservable gauge changes proportional to x'^ or require the absolute 
position of particles have been omitted. Such interactions include both space- and time-like four vector 
interactions, and under special relations between the a's and (3' s are called electromagnetic-like interactions 
due to how the potentials share similarities with the related time and space components of vector potentials 
in classical electromagnetism. Using the variables and , these can be incorporated into tt^ — pi — Ai 



A"^ 



ai(a;-L)Pi +^i(a;_L)P2 > 




a2{xi_)p2+ P2{xl^)Pl■ 



m 
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as 



= £2?" - G2P^ , 



(85) 



Thus 



TTl -pi 


= -^1 


= {Ei 


-ei)P- 




-l)p 


Al 


= (£1- 




-(Gi- 






7r2 -P2 


= -A2 


= {E2 


- S2)P - 


-(G2 


-l)p 


A2 


= (£2- 


E2)P 


+ {G2- 


l)p . 





When Ai is time-like (parallel to P), Gi — 1 = G2 — 1 = and 

Gi = G2 = Gtl — 1 • (87) 

One assumes that 

Gi = G2 = G = exp(g) (88) 

for general vector interaction. This implies that E\ , E2 ,and G depend on two invariant functions, which will 
be called A and V. A is redefined instead of G as a basic unit in the potential by choosing G = G{A). The 
most general choice for Ei is then Ei = Ei{A, V), so 



which modifies the free expressions 



Si(A,V)P^ + G(A)p^ , 

E2{A,V)P^' - G{A)p^' , (89) 



p1 = siPi'+p^, 

p^2 = £2P''-p^ (90) 
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the case of pure time like invariant V = and 



E,{A) =G{A){e^-A), (91) 

with 

corresponds to what will be called the electromagnetic-like vector interaction. With this choice the Hamil- 
tonian is then 

« -^2 = TT^ + mj = -El +ml + G^p^ - 2E2Gp ■ P « -Sf + ml + G^p^ 
= GHp^ - (ei - Af + ml/G^) = G^p" {s^ - Af + ml/G^) 

= G^{p^ + 2e^A-A^ -b^) (93) 
where Eas. ipT]) and have been used to obtain the result in the last line above. Thus one arrives at 

= (G' - l)(p' - b^) + G\2e^A-A^) (94) 

so that similar to Ea.(|55| 

H^G'^ij)'' +ml-{e^- Af). (95) 
Note that since £\~ E\ — e^ — E^ one still has the hyperbolic Eq.(|77l). but now Eq. (|9T1) implies that 

J = -Q. (96) 
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2.3.3 Scalar and Electromagnetic-like Combination 

If one combines both scalar- and electromagnetic-like interactions, one has 



Ui = ttI + MI = -EI + MI + G^p^ + 2EiGp-PK-El + Ml + G^p^ 

« -^2 = TT^ + Ml = -El + Ml + G'^p^ - 2E2GP ■ P w -El + Ml + G^p^ 

= G'^ip' -{ei- Af +Ml/G'')^G'^{p' -{e2- Af +MI/G''). (97) 

Following the analogy between Eas. (|83)) and (|95|) it is plausible, in light of Ea. (|60)) . to take 

n = G''{p' + {m^ + Sf-{e^-Af) 

= G^{p^-b^ + 2m^S + S^ + 2e^A-A^). (98) 



If one takes $^1 as given in Eg . ((94)) then 

$s = (G2-l)(p2_fo2)+G2(p2_&2^2m,„5 + 52 + 2£,„A-A2) 

= (2to„5 + 5*2) , (99) 

which, unlike Ea. (|59l) depends on the invariant A responsible for vector interactions. This equation and the 
last line of Eq.([97l) imply 

Ml = ml + G^ {2m^S + 5^) , 

Ml = ml + G'^ {2m^S + S^) , (100) 

and so 

Ml - Ml = ml - ml. (101) 

This again gives the hyperbolic solution of Eq. pS]) . We have derived the general form of the vector potential 
and then shown how it specializes to the time-like and electromagnetic-like forms for a one-body Hamiltonian 
and then continued to the derivation of the two-body Hamiltonian equations. In addition, we have then 
shown the combination of the scalar and electromagnetic-like vector interactions. 

What we have shown so far is entirely classical, so now that the potentials are completely defined in a 
classical formalism, we will quantize this system. 
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2.3.4 Quantum Constraint Dynamics for the Two-Body Relativistic System 

Now that wc have described all of the relevant interactions relating to the two-body Hamiltonians classically, 
it becomes pertinent to review Dirac's method of conversion to quantum mechanics from classical mechanics. 
First, the constraints derived above now become conditions on the wave function, i.e. "H « — > T-Lijj « 0. 
Secondly, Poisson brackets are now commutators. In the quantum case one must use 

p/i py 

K,p-]=i{v'--^^), (102) 

as a check to verify that the commutator ['Hi,'H2] (the quantum compatibility condition, equivalent to the 
classical {Wi, 'H2}) on the wave function vanishes. This ensures that the quantum versions of the constraints 



HiiJ = , (103) 
H2^|J = 0, 

are compatible wave function conditions, and 

P-}jwO-S^P-pV = 0. (104) 
Now we see how this applies to each of the interacting potentials. 

2.3.5 Quantum Constraint Dynamics for World Scalar Interaction 

This is a brief review of the methods for adapting the above section of the scalar interaction into the quantum 
realm. By setting <i?s = 2mwS + S'^(the scalar potential), the condition 



gives 



and the quantized Hamiltonian is 



P-p«0, (105) 



P-p«0-i>P-pV' = 0, (106) 



7^ = p2 + + $s « ^ (p2 + (m^ + Sf)^ = 0. (107) 
The eigenvalue form of the Hamiltonian above is, in the center of momentum system, 
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(p2 + 2m^5 + 52)V = 6^(w;)V. 



(108) 



Comparing this to the nonrelativistic energy expression 



(r) + 2mU (r) ip (r) = 2mEtp (r) , 



(109) 



one can see that 2mwS + and b'^{w) are analogous to 2mU (r) and 2mE, respectively. The scalar 
interaction has then been properly defined in the quantum realm, so we continue to a description of the 
vector interaction in a quantum formalism. 

2.3.6 Quantum Constraint Dynamics for World Vector Interactions 

To construct the quantum versions of the 'Hi's we must maintain the minimal substitution form of tt^ . 
Classically, the underlying scalars {S, V, A) may depend on x'^,x± •p,p^,and w"^ or x']_,P,p'^, and w"^ where 



is the invariant square of the relative angular momentum. A scale transformation on the wavefunction can 



so that has only time-like eigenvalues. Assuming that the Ei [= G{ei — Ai)] are Hermitian, then starting 
with the classical expression equation 



= (r X p)^, (in center of momentum system) 



(110) 



eliminate linear dependence on x± -p. The system variables p and P are well defined if the space is restricted 



n>l = Ei{A)P''+G{A)p^, 



(111) 



the quantum Hermitian version of which is 



K)h = 




(112) 



These equations combined with 
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E2{A) = G{e2-A2), 



(113) 



give us the Hermitian forms 



G 



< = G 



p''{e2-A)-'[^-^y^^g 



(114) 



where Q = \n.G. Using the quantum brackets 



(115) 



One verifies that the commutator [Hi , '^2] on the wavefunction vanishes and that the quantum analogs of 
the classical constraints {Hi = 0) become conditions on the wave function 



^iV = (tt? + m?)V = , 

H21P = {nl + mDip = . (116) 

Due to the two-body system being an isolated one and ijj being an eigenfunction of the total momentum, P'^ 
becomes a constant of the motion. The difference constraint 

(Wi-W2)V = 0, (117) 

then gives 

P-p^ = 0, (118) 

which is now a differential equation as p^^ = —id/dx^ in the coordinate representation. The quantum 
counterpart to the remaining independent constraint {H « 0) is Hip = for the system wave function. In 
the cm. system this yields 
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(119) 



When separate time and electromagnetic-like vectors and the world scalar interaction are used then our 
quantum constraint equation can be written as [4] 

G2[p2 - (e^ - Af + 2e^V -V^ + (m,, + Sf + -VG ■ p - ^-V^G - ^(Vg)^]^ = . (120) 

t 2 4 

If one lets (f> — Gip, the above equations become 



(£t, 



Af + -V^G 



= 0, 



(121) 



and[4] 



p2 - (e„ - A)' + 2£„V - + (m^ + Sf + ^V^G + \iVGf 



, 



(122) 



respectively [2S], which is the general Klein-Gordon form for two spinless particles. The terms with G can 
be regarded as quantum recoil corrections. That is, they would vanish if one of the particles becomes 
very massive. This completes the review of the relativistic quantum mechanical fundamentals of constraint 
dynamics for spinless particles, so we will now continue on to the relativistic forms for spin-one-half particles 
(i.e. the Dirac equation and how we arrive at the two-body Dirac equations). 



2.4 One Body Dirac Equation 

Since the potential framework that is used to extend the above wave equation for the two-body Klein-Gordon 
equations to the two-body Dirac equations in the constraint formalism has been assembled, it is helpful to 
briefly review the one-body Dirac equation in order to see how spin dependence comes into play. Defining 
the operator S aj3 

S=j,ij-p-Ep + M), (123) 



the Dirac equation can be written in the form 



= ■ p-E^ + M)i; (124) 



*The 75 here does not affect the physics but simphfies the compatibility conditions in the two-body case. 
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where E = e — V , and M = m + S, accounting for the scalar and vector interactions to be defined later. 
The 4x4 matrices are [2] 

I rr n \ 

(125) 
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:) 















With the 7g matrix 
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1 

1 0, 



and the four component spinor written in terms of two component subspinors 



^1 

^2, 



Then the Dirac equation becomes 



{-a ■ p)i;, + {E + M)i:^ = 0, 
(-(T.p)7/i2 + (S + M)Vi =0, 



(126) 



(127) 



(128) 



(129) 



(130) 



by "splitting" the matrix equation into two linear ones. If one solves the second equation for ?/'2 and 
substitutes into the first, one obtains 



[p2 -E^ + M'^ + i{(T-V In {E + M))cr • pj^/^i = 0, 
and similarly solving for t/'j from the first and substituting into the second 

[p2 -E'^ + -i{a-V In {E - M))ct ■ pjV'j = 0. 



(131) 



(132) 



If one expands using an identity of Pauli a matrices {<JiCfj = 5ij + i£ijk<^k) and take V ln(X) = ln'(X)f , 
one arrives at 
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[p2 -E^ + + i In' {E + M)y ■ p 



In' {E + M) 



r 



[p2 -E^ + -i In' {E - M)r ■ p+ 



In' {E - M) 



(133) 



r 



and thus we can see how the reduced Dirac formahsm reveals the famihar spin dependence with the — term 
and Darwin interactions, which together give the correct spectral results for all orbital angular momentum 
for Hydrogen[33]. This well-known result is explicitly displayed here to demonstrate in a simple form what 
the more complex two-body Dirac formalism does later. 

Now that the familiar formalism[14 for one-body Dirac equations has been reviewed, we can discuss 
two-body Dirac equations. 

2.5 Two-Body Dirac Equations 

This section presents a review of two-body Dirac equations j25l . The Dirac equations for two free particles 
can be written in the form 




(134) 



The theta matrices are defined as [35] 




(135) 



and these satisfy the anti-commutation and commutation relations 
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[^5^,05^]+ = -1, (136) 

[0,f,0^^]_ = 0, 2^ J, 
[05^,^5J]- = 0, l^J, 



In terms of variables defined earlier 

Pi=eiP + p, P2^e2P-p, (137) 
where P — P/w, Ea. p34p can now be written as 

Sioij = {ei ■ p + ei0i ■ P + miO^i)^ ^ , 

W = i-e2-p + e292- P + m2e52)^j = 0. (138) 
Using Eq. (|136p one can show that 

[5io,^2o]^ = 0. (139) 
The projections of these theta matrices satisfy 



[e,-p,e^,-p]+ = 1, 

[0,-P,C]+ = 0, (140) 

with 9^j^ = Oi^{ri^^ + P^P'^). These theta matrices are Dirac matrices modified in such a way as to ensure 
that the operators Siq and ^20 are the square root operators of the corresponding mass-shell operators, that 

is Sfo ^-UpI+ O and ^lo ='Wl + "^ifl 
The difference 

^In the ordinary Dirac matrix formalism, the Klein-Gordon operator (j>\ + m\) is obtained not from the direct square of 
Dirac operator 7i ■ pi -I- m\ but from multipUcations by the related operator — 7j ■ pi -|- m\ . 
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(141) 



leads to the equation 

P • = i [wiei -82)- {m\ - m^)]^ = 0, (142) 

which is the same constraint on the relative momentum that eliminates the relative energy in the center of 
momentum frame that appeared in the spinless case. 
Covariant Dirac a and /3 matrices are defined as 



/?. = -1,- P = 29^,9, -P, 

at = 2etJ,-P, (143) 

along with the E matrix 

= 2\/2i05.^»--P^±^, 1 = 1,2. (144) 

These take on the simpler forms — (0, OLi) and — (0, Si) in the center of momentum system for which 
P — (1, 0). In that case the free two-body Dirac equations take the more familiar forms of 

0^^{l^ ■ P - P,e, + m,)iP = ; i = 1, 2, (145) 

analogous to Eq. (|124p . or 

h^-V- P,e. + mi)i,^Q- 1^1,2. (146) 

If one makes the naive substitutions for scalar interactions defined by Eq. (1431) , the Dirac equations would 
become 



^iV' = (6ii-p + ei6ii-P + A//i6'5i)V' = 0, 

= (-02-p + e2e2-P + i\/2052)V' = O, (147) 

which are then not compatible, that is 
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[Si,S2]-tp = [01-p,M2952]- + [Ml051,-e2 -p]- 

= -i{dMi ■ 91652 + dA'h ■ 62951)1^ 7^ 0. (148) 

Using supersymmetry arguments and the third law condition of Ea. (|45p and Eq. (j46p . one can obtain the 
compatible spin-dependent constraints Si and 52 below[55]. 



Slip = {ei-p + ei9i-P + Mi95i-idC-02952O5i)ijj = O, (149) 
S2i> = {-92-p + e292-P + M2952+idC-9i95295i)^ = 0, 

which commute with each other. These equations are known as the external potential forms, due to how 
each particle's potential is created by the other particle. They satisfy 



[5i,52] = 0, (150) 

provided that 

ac^?Mi.?Mi. (151) 

M2 Ml ^ ^ 



The strong compatibility due to a supersymmetry is responsible for producing the spin-dependent recoil 
terms depending on dC. It can be shown that these terms vanish when one particle is infinitely massive 
relative to the other (Mi >> M2 or vice versa, giving dC = ^ = ^ = or dC = ^ = ^ = 0), 
thereby recovering the one-body Dirac equation. The operators Si and 5*2 satisfy the relation 

(5? - - -^{pI + ml -pI~ m^)^ = . = 0, (152) 

thus showing that the relative momentum remains orthogonal to the total momentum. 

The above difhculty with compatibility from Ea. (|148|) can also be overcome if the scalar interaction is 
replaced by vector, pseudoscalar, pseudovector, or tensor potentials as described below. We first rewrite 
the Dirac equations Ea. (ll49p for scalar interactions in the hyperbolic form of|2H] 
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Si^ = (cosh(A) Si + sinh(A) 82)^- = 0, 

S2tp = (cosh(A) S2 + sinh(A) 81)^- = 0, (153) 



where 



and 



Siilj = (5io cosh(A) + 520 sinh(A) )V' = 0, (154) 
S2V = {S20 cosh(A) + S-io sinh(A) )^ = 0, 

A = -95i052C{x±). (155) 



Here, 81 and 82 are auxiliary constraint operators. 
It turns out that for all arbitrary A satisfying 



[P-p,A]V = 0, 

Aix) = A{x^), 
P-pi) = 0, (156) 



this leads to weak compatibility on both sets of constraints 



[Si,S2\i> = 0, 

[81,82]^ = 0. (157) 



The function A is defined as 



for scalar interactions, 



Ac = -Ce^ie^2 = = -751752, (158) 



Aj = JP-e^P -92 = 02^ = fi^p2^0^, (159) 
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for time-like vector interactions, and 



(160) 



for space-like vector interactions. Here, C^J^, and Q are all invariant arbitrary functions of xj^ related to 
each interaction. For the axial interactions, see [28] . |35j . 

2.5.1 Quantum Analogy and Reduction to Schrodinger-like Form 

Similar to the one-body case, this section discusses an analogy of the two-body for converting to Schrodinger- 
like equation. This reveals one of the advantages to our methods in that methods for solving Schrodinger- 
like equations are very well known and so if we can convert any formalism into this form, a final solution 
becomes much more feasible. For the combined scalar, time-like vector, and space-like vector, the equations 
Eq. (|154ll53p take the form below after bringing the free Dirac operators through to the right p5].[55] 

5iV^ = (G01 • p + EiOi ■ P + Afi05i + i^{e2 ■ dgOi + 02 ■ dJ02 - ^2 • dCOi))^: = 0, (161) 



S2ip - {-GB2 ■ p + E2B2 ■ P + 1/2^52 - i^{Oi ■ dg03 + 01 ■ dJ02 - 01 ■ dCOi))^j = 0. (162) 

In the two-body case, the wavefunction is now a 16 component spinor, with each individual having four 
components 

■02 



(163) 



The end result of the matrix multiplications of Sitjj = is a set of eight simultaneous equations for the 
Dirac spinors V'li "02 7 V'si ■'/'4- There is a set of four simultaneous equations from the constraint of the two- 
body Dirac equations that Sitjj = and another set of four simultaneous equations from the second of the 
constraint two-body Dirac equations 324^ = 0. In order to fully solve these equations, we define the terms 



4'± = "01 ± '04> 
X± = '02±V'3: 



(164) 
(165) 



D++ = e^{ai-p - ^-02 ■d{J-C + gai ■ ^2)), 



(166) 
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D++ = e^{cj2-p - '-cJi ■d{J-C + Gai ■ aa)), (167) 

=e^{ai-p- ^(72 ■diJ + C + Gai ■ tJa)), (168) 

D-- = e^{<Ji .p+^a2-d{J + C- Qai ■ cia)), (169) 

El = £1 cosh(j7) + £2 smh{J), (170) 

E2^e2Cosh{J)+eismh{J), (171) 

Ml = mi cosh(/:) + TO2 sinh(/:), (172) 

AI2 — 1712 cosh(£) + mi sinh(£), (173) 

and work out the matrices to obtain 

=i?iX+-Mix_, (174) 

D+-^+^E2X+ + M2X-, (175) 

D^+X+=Ei^+~Mi^_, (176) 

D^-X- = -Ei(j)_ + A/i0+. (177) 

These equations can be solved through a system of substitution and ehmination similar to what was done 
for the one-body Dirac equation (Ea. (ll34p '). since there are now four equations and four unknowns. The 
result is 1551 



(i?2-Mf)0+. 
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Further expansion of terms reveals [3 5) 



[p2 + 2m^S + + 2e^A - 

^ E2Mi+EiM2^ ^ 2 ^ 4^^-' ' E2Mi+EiM22^ ^ ^' 

i>j_+_CT2) . , _ 1 £^2^2 + EiMi _ _ L-{ai- C72) E2M2 - E^Mi , _ , 
r 2E2Mi+EiM2^ 2r E2M1 + E^M2^ ' 



^ ^4 ^ 2r 2^^' 2^^ ^' E2M1+E1M2' 
+ai.ra2-r(-V£--£ +gL - -C {C ^ G) ^ ) 

.■i,r,r\U - ( - \,^ir £^iA^2 - £^2^1 i • (^i x fJa), , 
+ + G) (^1 ■r(T2-p + (J2- rai ■ p) + -(£ - G) ]0+ 

= + 

= &'H'/'+, (179) 

We work in the cm. frame in which P = (1,0) and f — (0,f). The final four component wave functions 
V'i , ?7± that appear in Eq. (|179p are defined by [5] 

(/)_!_ = exp(J' + /Ccri-fcr2-f)'(/;_|_ = (exp J^)(cosh/C + sinh/Ccri-f(T2-f)'(/'_|_, 

X± — exp{J^ + ICcri-rcr2-r)rij. — (exp J^) (cosh /C + sinh/Ccri -f 0-2 -f)??.!., (180) 



in which 



= 77 log 



2 £2^71-1 + ei"i2 
V=E2Mi+EiM2, 

jC^i^—^. (181) 



In analogy to what occurs in the decoupled form of the Schrodinger equation for the individual single particle 
wave function, this substitution has the convenient property that in the resultant bound state equation, the 
coefficients of the first order relative momentum terms vanish and we are left with a more concise form 

(p^ + $12 - fo^)?/' ={p^ + (m„ + Sf - (e„, - Af + ^D + ^- {(Ti+(T2)^so + (Trf(T2-fL • {rTi+(T2)^sOT 

+ (Tl-(T2^SS + (3(Ti-f(T2 •f-(Tl-(T2)$T + L • {(Tl-(T2)'t gQD + iL-CTi X Cr2<l>SOx}V'+ = 0. 

(182) 
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Here, J ^ Q and £ arc arbitrary functions of related to each interaction, as defined in Eq. (|158p . Ea. (|159p . 
and Eq. (jl60p and Q = — J . Thus is derived a two-body Schrodinger-hke equation for scalar, time- like and 
space-like vector interactions. The 2m^S + 5*^ + 2e^A — term is the classical interaction potential 
term (which also appears in the spinless Klein-Gordon equations), the L ■ (g\ ± a^) terms represent effective 
magnetic field dipole moments and Thomas precession, and a\ ■ rai ■ f, o\ ■ and Oi ■ faj ■ p terms arise 
from dipole-dipole interactions and their relativistic corrections. The main focus of this work has been to 
derive a similar equation for the three-body baryon system. 

This is the framework for the two-body system in a fully relativistic formalism, so from here we go to 
larger systems. Sazdiian|ll) has done considerable work on the A^-body system, which will be reviewed 
shortly. Although he does not deal with spin dependence with as much detail as done here, he does provide 
a very useful framework for the A^-body problem in a constraint formalism. 

2.6 Two Body Dirac Equations: Explicit Forms of the Potentials 

Since the forms of the potentials in the three-body case are similar to those in the two-body case, it is of use 
to (very) briefly describe the two-body interacting potentials and how they affect the wavefunction. This 
section then contains a review of how the operators of the tensor, spin-spin, spin-orbit, spin-orbit difference 
and spin-orbit exchange work on a {jlsn\ state coupling, where 71 is a radial quantum number. In the 
three-body case there will be two n's, one for each relative coordinate, but here there is just one. 
First, our quasipotential (energy dependent effective potential) is defined by 



^ = ^SI + $D + L ■ (cri+(T2)$S0 + (Ti-fcr2-fL • (cri+(T2)$50T (183) 

+ cri-(T2^ss + (3cri-fcr2 -f-cri •cr2)$T + L • {cri-(T2)^ ^OD + jL-o-iXcr2$sox, 



where the potential term^ 

^D, ^SO, ^SOT, ^SS, ^T, ^SOD,^SOX 



are all collections of two-body terms depending on the masses, distances between the two particles, energies 
of the two particles and the energy of the system. The explicit forms of these are therefore not important 

®The subscript on quasipotential refers to Darwin. It consist of what are called Darwin terms, those that are the 
two-body analogue of terms that accompany the spin-orbit term in the one-body Pauli reduction of the ordinary one-body 
Dirac equation, and ones related by canonical transformations to Darwin interactions [7, momentum dependent terms arising 
from retardation effects. The subscripts on the other quasipotentials refer respectively to SO (spin-orbit), SOD (spin-orbit 
difference), SOX (spin-orbit cross terms), SS (spin-spin), T (tensor), SOT (spin-orbit-tensor) 
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to the current discussion of the operators and so will be left in this form for simplicity's sake. The spin 
independent and Darwin terms have no spin operators and so when used on a {jlsn\ state they just give 



{jlsn\<S>si\jl's'n') ^ Sw5ss'{n\<i>si\n'), 
{jlsn\^D\jl's'n') = SwSss'{n\^si\n'). 

These are the only simple terms though, as the spin-orbit gives 

{jlsn\L ■ (cri+cr2)$so|i^'sV) - + 1) - l{l + 1) - 2]Swd,,,S,i{n\<i>so\n') . 

As we will show later, while in the two-body case this Z,j,and s are for the entire system, in the three-body 
problem it is just for each pair of particles and so this spin-orbit function requires additional (and extensive) 
manipulation in order to reach a completely coupled \JLS) state for each set of particles. The emphasis 
here is important as this is the main difficulty in going from the two-body formalism to a three-body one, 
in this work as well as others. 

The tensor and spin-orbit tensor terms allow for coupling of different / states as well as identical / states, 
as shown ( ^49_ ) 

Tensor: 



2? + 4 

Oj -f lln|(3(Ti-f(T2 •f-cri-cr2)$T|jj + 11?^') = -— — -(n|$T|n') 

2j + 1 



{j3 + Un\{3(Tvr(T2-r-(Ti-(T2)^T\jj-nn') = 6 ^^^ {n\^T\n'), 



2j + 1 

0j-ll7i|(3cri.f(T2 •f-o-i-cr2)$T|ij-lln') - -|^(n|$T|n'). (184) 



Spin-orbit Tensor: 



{jj + lln|cri-fcr2-fL • (cri+<T2)$soT|jj + Hn') = -— — [-2j - 4] (n|$50T|n'), 

2 J + 1 

{jj + Un\avr(T2-rL ■ {(Ti+a2)^soT\jJ - Un') = [-2j - 4]2^^^^^^{n\<^soT\n') , 

2j + l 



{jj - lln|(Ti.fcr2-fL • ((Ti-f^2)$soTbj + lln') = [2j - 2]2 ^ ^^^\ n\^ sotW) , 



(jj-lln|crrf<T2-fL-(cri+cr2)$sOTbj-lln') = —^[2j-2\{n\<^soT\n'). (185) 
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2i + i 

2j + V 



The spin-orbit difference and spin-orbit exchange only allow couplings between different spin states 



071n|L.(<Ti-^2)*5O^|j7'ln') = 



071n|L.K-^2)*soDlj''0^0 



[{1 - (j - + m + I + 2mn\^soD\n') 



0m|L.K-^2)*soDlj''l^0 



[{1 - (j - 0'}{(J + 00- + I + 2)}](n|*soD|n') 



070n|L.K-^2)*soDb'''0^'> = 



(186) 



{jlln\iL-(TiX(T2'^sox\jl'ln') = 



[{1 - 0' - 0'}{(J + 00' + I + 2mn\^soD\n') 



{jlQn\il.-a-i x cr2$sox \jl'ln') 



[{1 - (j - Ifm + m + I + 2mn\^soD\n') 



(j70n|iL-criX(T2$soxb7'0ri') = 



(187) 



The exact derivations of these potential terms are done in ([H HJ |3]). Since one of the goals of this work is 
to compare essentially the same methods that worked well for the meson spectrum to the baryon spectrum, 
we use these same potential terms in mostly the same form as they appear in the pure two-body case. The 
two-body operators are therefore loosely defined and described in preparation for their adaptation to the 
three-body potential. Now we will derive definitions for the scalar and vector potentials used in our model 
and from that define the two-body potentials ^Dt'^so, ^sot, ^ss, ^t, ^sod, and ^sox- 



The authors of [52] have used a sophisticated form of the static quark potential developed by Adler 
and Piran [3^, one that has ties at all length scales to field theoretic data and from this obtained good 
agreement with the quarkonium spectrum from experimental data. However, it is much more common in 
nonrelativistic treatments to use the static quark Cornell potential ^4T for potential model studies. 



as in 133] . Although not displaying asymptotic freedom, it does give the dominant Coulomb- like behavior 
as well as the linear quark confinement. Early on a model was proposed by Richardson for a static potential 
which both depends only a single scale size A and interpolates in a simple way between asymptotic freedom 



2.7 Explicit Forms of the QCD Model Potentials 



V{r)^-—+br, 



(188) 
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and linear confinement '44'. Richardson's model for the static interquark potential in momentum space is 



arising from the assumption that 

3q 



n^)--^^' (190) 



(including the color factor —4/3). It is important to note that this is for a qq singlet state for the meson 
spectrum. In order to properly account for asymptotic freedom, we must have q^/A^ >> 1, which gives 

On the other hand, the property of linear confinement requires that for Ar >> 1, V{r) cx r or equivalently 
that for q^/A^ << 1 one must impose ^^(q'^) ^ q ^ . The interpolation of Eq. (|189l) is not tied at all in the 
intermediate region and only roughly tied in the large r region to any field theoretic data. Nevertheless it 
provides a convenient one-parameter form for the static quark potential. In coordinate space it has the form 

n.) = ^-^^, (192) 

^ 27 27r ' ^ ^ 



rrrQ 



where /(Ar) is given by a complicated integral transforrrlJ that displays the asymptotic freedom behavior 
for r — > of 

/(Ar) ^ --i-, (193) 
In Ar 

while for r — > oo, 

/(Ar) ^ 1. (194) 



A simpler model for the potential function /(r), which we apply in this model and one which displays 
the same large and small r behavior ijf| 

_ 8^A^r 167r 

^"^^ 27 27rln(e2 + l/(Ar)2)- ^ 

It amounts to replacing Richardson's /(Ar) by 2/ln(e2 + l/(Ar)2), having the same limits. The forms of 



^In addition to the spin independent nonrelativistic model presented in 1441 see also a relativistic extension of it given in 

m- 

*An earlier coordinate space form that displays asymptotic freedom as well as linear quark confinement proposed in | 46l is 
V = (87r/27)(l - Ar)2/(rlnAr). 
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the scalar and vector invariant potentials 

87rA2 



S 



27 



A = '-^ ^ + — (196) 

27rlog(iv:e2 + ^) 47rr ^ 

are used to construct all of the individual $ terms. The explicit forms of these potentials that are derived 

from the vector and scalar potentials arc given in Appendix A. In the case of the baryons these are slightly 
changed due to a different color factor (— becomes — due to this being quark-quark and not quark- 
antiquark) to 

J — 



27 



Stt 6162 



27rlog(Xe2 + ^) Anr ^ ' 

and also of course there is no longer just one interaction but three, so r becomes ri2,ri3,or r23, depending 
on which potential we are currently discussing. 

The technique that Crater et al. used in the two body problem for finding the eigenvalues is called the 
Inverse Power Method and its application depends on the variables being separable. Unlike the two-body 

problem, the variables are not separable in the three-body problem. This requires the use of the variational 
principle, which in turn requires a basis. Therefore, now that the potentials are defined, we turn to defining 
a basis for the two-body system and from there extending it to the three-body. 



2.8 Gaussian Basis Functions 

Here we will briefly detail how our wavefunctions are used to construct the basis for use with the 
variational theorem. In the case of a nonrelativistic, energy independent Hamiltonian the standard way of 
using the variational principle is having the wavefunction vary by some alpha, i.e. 



{^{a)\H\^{a)) = E{a), 

— - 

da ' 

E{ao) > Eactual- (198) 
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In order to expand the basis into a more general matrix eigenvalue equation, we define. 



I*) = ^c„|vI/„), 

71 

{mm = ^c„c;.(vi.„.|i7|vi/„), 

n.m 

= ^c„c:;,(^'|*)=^c„c:;,B„™. (199) 



Since the basis we will choose for the baryons is not orthogonal, B is not the usual Kronecker delta. Using 
the method of Lagrange multipliers (subject to the constraint (^'1^') = 1) gives that 

^ ((^-li/l*) = ^c„(*,„|i7|^'„)-S^c„B^„ = 0, (200) 



dc] 

we arrive at the eigenvalue equation 

or in matrix form (where H and B are matrices, c is a vector and is a scalar eigenvalue) 

He = EMc. (202) 

In the non-relativistic two-body problem with tensor coupling, the s and d states {I — j + \ and I — j — 1) 
are mixed, so we need a mixed wavefunction 

1*) = ^C„+|*„+)+E^"-l*"-)' 

n n 

+ ^ l=J + l- (203) 

Using this * in Ea. (jl99l) gives 

n.m 

(*l*> = Y.^c*+c++c*^c++c*+c-+c*^c-), (204) 
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and as before, we take the partial derivative with respect to the complex conjugate of one of the c parameters 



dC* 



i^\H\^) - E(^\^) = H++C+ + H+-C- =EB++ct 



H — + H^^c^ 

mn ' mn^n 



= EB- 



(205) 



giving us 




= E 




(206) 



Now we must define a wavefunction for our basis. 

2.8.1 Gaussian wavefunctions and Infinite Interval Discretization 

This section describes how our wavefunctions comes about for our basis. The potentials in our Hamilto- 
nian have both short distance and long distance effects, so we need a basis wavefunction that can accurately 
account for that. We define a wavefunction in terms of some parameter a that determines the effect of the 
Gaussian wavefunction for short and long interactions. We then split the wavefunction into those two parts 
(short and long) and discretize it to a certain N value, from where we get our basis. The wavefunction is 
originally defined in an infinite vector space, so we must truncate it in order to work with it. 

Boris Kupershmidt, a mathematician. [50] has suggested a Laplace transform/Gaussian basis 



^(x) 



daq{a)] 



I , ax^ , 



(207) 



where t/i is essentially the Fourier transform of some function d. In order to work with this function, we split 
the integral into two pieces, one with boundaries from zero to one and and the other with boundaries from 
one to infinity, so that 



^(x) = ^ da9(a)^iyjexp(-|^) 

r . . N / 1 , ax\ 

+ A da'Z(a)y^V^^exp(- — ). 



(208) 



By replacing a with in the first half of the equation (so that the integral from to 1 now becomes 1 to 
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infinity) we get 



exp(- 



daq{a)] 



■ exp(- 



2a2 



and from there, replacing integrals with sums over arbitrarily large N, this discretizes to 



rt=l 



1 Tjx^ 



1 / 1 



1 / 8 



N 



n=2 ' n=2 



a"" V TT"' 



So, for A'' = 1 we have 



V'W = eJ^y^exp(-ii,), 



For N = 2 we have 



Note that the original wavefunction from N = 1 remains as the first term. This is true for all A''. 
For N>3 



n=2 



2 

2na?' 



1 /n3 



nx 



exp(-Tr:T) 



n=2 
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or more symmetrically 



V'(x) 

fn 
In 



2N-1 
1 



V 



l<n< N: 



n+l-N: N + l<n<2N-l 



(214) 



As we can see from the N = 2 case, the order of the matrix increases as 2N — 1. Each matrix element 
of the Hamiltonian matrix is constructed from the expectation value of the Hamiltonian with two of these 
wavef unctions. For example, for the N = 2 case, our general wavefunction |i/'„(/n)) is 



IV^i(i)), 



= 2^1^2(3))' 

= 3^|V^3(2)), 



(215) 



and thus we have the 3x3 matrix 



(216) 



As can be inferred from the values of /„ for n > 1, this basis allows the wavefunction to account for both 
the long-rage and short-range interactions of the Hamiltonian. Smaller /„ values-such as for n = 2 in the 
above example-allow for long-range interactions while larger /„ values (like the n = 3 wavefunction) account 
for the short-range interactions. 

In a similar manner, we can now also write our B matrix from Eq. (j202l) as 



2N-1 



V'(x) 



B„ 



E '^"V'nlx); 



d^xV':(x)V'„(x) 



{fn + fmW 

2a2 



) 



/.3/2 „3/2 
V Jn Ji 



Q f 3/2 /•3/2 



{fn + fmfl^ V {fn + fm? 



(217) 



44 



Thus we get an analytical form for the B matrix that remains the same regardless of coordinate transfor- 
mations. Note also that this becomes one in the case of /„ = /„ = 1, which is expected. This completes 
our review of the two-body formalism. Since we are attempting to reach a convergence point with as few 
Gaussians as possible, we do not necessarily include as many wavefunctions as is possible. So for N = 2, 
we only begin with two wavefunctions for each coordinate (giving a 4x4 matrix) and then go to three wave- 
functions (going to 9x9). Similarly, = 3 can have up to 5 wavefunctions per coordinate, but we only add 
one at a time in order to more quickly converge the energy eigenvalues. 

3 The Three Body Problem 

Now that we have completed our review of constraint dynamics for the two-body problem, we move on 
to the three-body one. Doing so is not a trivial task, as many pieces of our Hamiltonian become considerably 
more complicated, or at the very least, considerably different, from those used in the two-body problem. We 
will begin our review of the three body problem with a simple non-relativistic Hamiltonian with a harmonic 
oscillator potential. As is common for three body approaches, we set up a relative coordinate system from 
the coordinates of the quarks themselves. Prom there, we change the potential to a power law and add in an 
ad-hoc spin dependent term to acquire crude spectral results for our potential, the purpose of which being 
to familiarize ourselves and the reader with the larger problem at hand, as while this is a vastly simpler 
problem, it still has the same form as the fully relativistic, spin-dependent Hamiltonian. The actual spectral 
results are found using a minimization routine. 

After the non-relativistic results, we move on to Sazdjian's relativistic N body formalism that we then 
adapt to a three-body formalism and from there further modify to a Schrodinger-like form, as we did in the 
two-body. We then once again set up our relative coordinate system, but use a more direct and advanced 
method to do so in order to avoid mathematical difficulties for the intrinsically more complicated system 
(i.e. dealing with an error function as we do below is not difficult, but it would become a much bigger issue 
with a more complicated Hamiltonian). Finally, we modify the two-body potentials as needed in order 
to fit a three-body problem. Some of these potentials are simple and require little other than changing 
an interaction between two particles to three interactions between three sets of two particles, while others 
require a nearly complete overhaul of the original potentials from the two-body system. 

3.1 Nonrelativistic Three Body Problem 

Before going to Sazdjian's work on the N-body relativistic approach (and from there to our three-body 
case), we begin with a simple non-relativistic three-body approach. The purpose of this is to provide a 
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much simpler analogue to the work done in the following sections in order to more readily facilitate an 
understanding of the workings of the system. Therefore, we first solve the non-relativistic problem for a 
three-body system using relative coordinates for a simple harmonic oscillator. Then, we replace the oscillator 
potential with a power law potential and add in an ad-hoc spin-dependent term and numerically solve the 
equations. We limit our work here to the octet and use a chi-squared minimization routine to acquire the 
best fit for the baryon octet. The non-relativistic Hamiltonian is 



H 



P? 



p1 



P3 



+ V{vi,r2) + V{r2, ra) + ^^(ra, ri). 



2mi 2ni2 

which in the simplified case of three coupled harmonic oscillators becomes 



(218) 



H-- 



9 9 

Pi , P2 



^ -r-^+iA3(ri-r2)' + iA3(ri-r3)' + iA3(r2-r3)' 

2mi 2m2 2m3 2 ^ ' 2 ^ ^' 2 ^ ^' 



(219) 



where is the usual spring constant with A having the dimensions of energy. Using the canonical trans- 
formation [9] 



P 
X 

R 



r2-r3, Pp = P2-P3 
Mm2 



(m2 + mz)mi 
miTi + m2r2 



r2 



(m2 -I- 1113)7711 
7713V3 



ra, Pa 



m3P2 - m2P3 
m2 + ma 



TOl + 7712 + 7773 



PI+P2+P3 



(220) 



where 



we arrive at 



M = 7771 -I- 7772 -I- 77I3 , 777^ 



7771(7772 -I- 777,3) 
(mi -I- 7772 + ^,3) 



777 2 7773 
7772 -I- 7773 ■ 



(221) 



H 



2M ^ 2777 A 



_ 

2mp ~^ 2 



(7772 -I- 7773)2 



7772 + 7773 



(7772 - mz)p ■ A+2A^}. (222) 



As a side note, similar to in Isgur's work, in the simplified case of 7772 = 7773 (most of the octet has two 
masses equal, so this is reasonable for this case), the Hamiltonian becomes 



H 



Pi 



Pi 



2M 2777A 2777p 4 



^A%p'+2y 



(223) 



It is interesting to note that this substitution bears a resemblance to a quark-diquark system in that in 
setting up our coordinate systems two quarks are favored over the third, but the similarities end there as 
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everything past that point treats the system as symmetrically as possible. Also, the oscillator case is similar 
to a typical phonon problem; however, since this case is just a simple example of what will be a general 
power law potential, it is not particularly relevant to the discussion at hand. 

This Hamiltonian Ea. (|223p is fairly simple to solve (for the ground state) using the Gaussian wave function 

3/2 3/2 

A) = e-a>V2e-o^.AV2^ (224) 

where 

ap = (3A3mp)i/4, (225) 
ax = (3A3m,)i/4. 

Now we will begin to build upon the simplest state by adding in spin dependent terms. The spin terms 
typically have the form 

V^— s,.s,V2v(r„r^.), (226) 

i<.j ■> 

(with c an adjustable parameter) that modifies the spin independent two-body potential and will be treated 
as three two-particle interactions. This spin term should not be taken as a derived term, but it is important 
to note that it bears some similarities to one that is derived later for Two-Body Dirac equations. The 
purpose of adding this spin term is simply to see how it might be treated in this much simpler case than 
the one with which we will actually be working. The results of working out the quark model expectation 
value Si ■ Sj terms are shown in Appendix C, along with the quark state vector structure for the baryons 
(Appendix B). It should be noted that though this addition does technically make the problem relativistic, 
we are still using the non-relativistic kinetic and potential terms from earlier, so for all intents and purposes 
it is a "non-relativistic problem with spin." Moving from the oscillator to the general case, the Hamiltonian 
now is 



Imi 21712 2m3 



+ 5]^s,.s,V2v(r„r^-). (227) 



i<,j 



'.mirrij 



For the purpose of further exemplifying the techniques we will use for the relativistic case, we choose V = 
iA(A|ri— Tj 1)^, which with the substitution from Eq. ()220p makes the Hamiltonian 
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2M 2mp 2mx 2 ' m2+m3 ' ' ma + ms ' 

+ Y.^.^^-^^^'^^''^^'^3)^ (228) 

where p is an arbitrary general power. The kinetic term can always be worked out exactly and is, working 
in the center of momentum system, 

(^|T|*) = ^ L-<p\-^'^>^" ^e-<P' (229) 

2mA 



as can the first potential term of |\/2p|P, yielding, for integer values of p, 



r 13 5 p+1 
2 2 2 ^ 2 

The other two potential terms do not work out analytically in the general case though and require some 
numerical work. Through the substitution 

a' = — — p-A, (231) 

m2 + ms 

one can move some of the more difficult terms (the I — p — Al^" and I A ^ — pp terms) into the 

Gaussian wavefunction, which gives us 

(*l|-;;rr;^'^-^l'l*>J^'^' = J^'^' / exp(-ay)exp(-ai(--^p-A')^)AVpd3A'. (232) 
The term in the second exponential above results in 

- - ^^P-AT = -«i(A'^ + p^+^^Apcosg), (233) 

rtii + TO3 (m2 + mzY m2 + ms 

which will cause the Hamiltonian to integrate to an error function in A . This can be approximated and then 
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Table 1: Nonrelativistic model parameter values and baryon masses 



PcLi'ticlc 


Mass(MeV) 


]V[a,ss(]V[eV) (GxpcriniGntal) 




u 


317 




d 


322 




s 


495 






971 


938 


Constants 


value 


n 


973 


939 


A 


175 MeV 


s+ 


1146 


1189 


c 


1.22 (dimensionlcss) 




1148 


1192 


P 


2.05 




1150 


1197 




^0 


1323 


1314 




1325 


1321 


A 


1125 


1115 



the remaining integral over p can be done numerically. The other term (|A — p|^)is done in a similar 

fashion and also results in an error function. We use a Monte Carlo modeling approach in a chi-squared 
minimization routine. The numerical results are as follows for the baryon octet 

So, we have now completed a simplified version of the three-body problem and obtained some crude 
spectral results. We have used a relative coordinate system for an oscillator basis at first and then changed 
our potential to a power law with spin-dependent terms. From there, we numerically solved the problem and 
used a minimization routine to acquire a best fit for the parameters (the quark masses and the constant 
A). 

3.2 Sazdjian's N body formalism and the Three Body Problem 

This section will focus on reviewing Sazdjian's work on the two-body and N-body systems[TT]. We 
describe his derivations for the A^-body problem and distill them down to a three-body formalism that we 
can then use for our specific problem. By using similar constraints to those defined in the two-body case 
for compatibility, he derives the same eigenvalue equation as in Eg. ([551 

{pI± +^a- b^MXa±) = [Ha - 6V](^a±) = 0, (234) 

where pj_ ^ p — PP ■ p is the transverse momentum. 

Sazdjian applies the covariant formalism with constraints to the N-particle case. The compatibility 
condition is then 

[Ha,H;,]* = (a,&= 1,...,7V), (235) 
which are N(N — l)/2 in number and give conditions on the interaction potentials ($a)- Later it will 
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be necessary to specialize Sazdjian's work to the three-body system so that it more directly applies to the 
three-body problem at hand. 

Sazdjian now includes all two-body interactions into each in order to construct separable interactions, 

i.e., 

$a=^$a6 + VFa, (236) 
6 

where Wa are remainders (in other words, three-body potentials) to be specified later. These potentials are 
symmetric {^ab = ^ba) and are solutions of the compatibility equation 

\pl + ml + ^ab,pl+ml + Ha]=0, (237) 



which since = can be written as 



\pl-pl,<^ab]= 0. (238) 



The two-body potentials $ab cannot alone satisfy the global compatibility condition. Therefore Sazdjian 
adds the Wa terms. The equations that must be satisfied by the Wa potentials are 

[plWa]-[pl,Wb] = -[Y^<^ac + Wa,Y,'^bc + Wb] {a,b = 1, N). (239) 

a, 6 a, 6 

Sazdjian makes a few remarks at this point. First, the total potentials would not have a closed form like 
in the two-body case (where $ is a function of x_l). Second, in the quantized two-body case, the potentials 
would become non-local operators, due to the two-body momentum operators Pab = Pa + Pb no longer 
representing the total momentum of the system and therefore not possessing corresponding eigenvalues. On 
the other hand the total momentum in two-body case is a number (an eigenvalue). He then abandons this 
approach and instead takes a simpler one that works only for confined systems where questions of correct 
cluster decompositions do not have to be addressed. 

Sazdjian finds that in working with the free N-body system, one can be guided by the simplifying features 
of the two-body system. Since the system must reduce to that of the free case in the absence of interactions, 
he found it useful to begin with the N-body system without any interacting potentials. In the free case the 
wave equations are 

ipl + ml)^{xi,..XN) = (a = 1, N). (240) 
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Let us introduce longitudinal and transverse momenta, which are defined by 



Pa = P\\a+P±a, 
P±a = Pa+Pa-PP, (241) 



where 



P = Pl+P2 + ---Pn 

and the cm energies of the individual particles which are 

P\\a = -Pa - PP = SaP, 

Pfla = -4- (242) 
By subtracting the two wave equations from each other we obtain 

{pl-pl)^ = {ml-ml)^, (243) 



or 

(pf|(, - pfi J* = {el - el)^ = [{ml - ml) - {pl^ - p^^)]* {a, b = 1,2, 3). (244) 

Even in the free case, however, the longitudinal momenta (p||) cannot be taken as the analogs of the lon- 
gitudinal momenta for the two-body case since plj_ — plj_ ^ 0. Therefore, one must define generalized 
longitudinal momenta. 

Sazdjian defines U as the operator of the canonical transformation 

U = exp[i^XaLkaL], (245) 
a 

where kaL are functions of the transverse momenta, total mass squared, and free mass, satisfying the condition 

Y,kaL= 0, (246) 

a 
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and XaL is a function of the particle position 



XaL =Xa-P. (247) 

The operator U transforms the Si into e^, where 

el = Ea + fcaL, (248) 

The operators ki are chosen such that the new (transformed) wavefunction is an eigenfunction of with 
eigenvalues that satisfy ( dropping the primes) 

el - el ^ ml -ml (a, 6 = 1, iV). (249) 

By dividing the above relation by (eq + £(,) and summing over all relations with respect to b for fixed a, 
Sazdjian arrives at N equations: 

Nea - E ^r'li^ = (« = 1: (250) 

where w is the total cm. energy 

w^Y.^b- (251) 

b 

These equations cannot be solved simply except in the two-body case, which gives 

{ml - ml) 
2ei - ■ = w, 

(ei + £2) 

2e.-^^f-^ = (252) 
(£2 + £1) 

which reduce down to Eq. ([29|) 

er-S2 = ^^^^, (253) 
w 

as expected. He shows it is possible, however, to find an approximate solution by using successive iterations 
in the general case. 

In order that Eq. (|249p be satisfied, the operators ka must satisfy the relations 

{Sa + kaL? -{eb + kf^L? = {ml-ml) + {pl^~pl^) (a, 6 = l...iV). (254) 
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By using a similar procedure to that for Ea. (|250p . he derives the relation 

{Pb± - Pa±) + (Sa - £b){kaL + hh) 



NkaL + y'"" ' ^"'7°" ' =0 {a = l...N). (255) 



b 



{Sa +£b + kaL + fcfoL 



An approximate solution for the k's can be obtained by keeping those terms which contribute in the 
non-relativistic limit. In that limit, Sa and pa± behave as c^, which means that ka behaves as c°. So, by 
using those approximations one comes up with the end result 



N 

|t[l + (^-Af)/(2m,m,Ecli2;b)] 



jV + jV^ h , ..n/.o^ ^.^N r^(« = l,-,^^)- (256) 



This equation gives the constituent cm. energies of each particle in terms of the total cm. energy and the 
masses. In the N = 2 case, we recover the exact two-body equations for the epsilons of 

w 1 (mi — 7712) m\ — + 

£1 — ' — 



£2 



2 2 [1 + (w — mi — m2)/(TOi + 7712)] 2w ' 

w 1 (m2 — mi) —m1 + m2+'w'^ 



2 2 [1 + (w — mi — m2)/(mi + m2)] ' 2w ' 

ml -ml 

£1-62 = , (257) 



and so our full N-body equation can now read as 



+ +m^]V^ = 0. (258) 

a=i Eb=i V2eb 

The equation (I258P determines the total cm. energy w in terms of the masses. Again, for N = 2, this 
becomes 

£1 + £2 £1 + £2 

which, since 



22 2 2 7,2 

£1 - mi £2 - m,2 = b , 

£l+£2 = w, (260) 

is the same as the equation derived in the two-body system earlier 

[^] + [^]^ = is2 + £i)^^ = (pI b')^P = 0. (261) 

WW w 



53 



Finally, he also derives the system of N Klein-Gordon equations 



N 



6=1 ^ c=l 



(262) 



Sazdjian finds that in the interacting case, the N non-independent wave equations are kinematic in nature 
and should not be modified by the interactions. For example, particle a "feels" an interaction potential $a 
that enters additively into its kinetic energy term by the relation 



Pla Pla + *a. 



(263) 



Ea. (|258p then becomes 



which in the two-body case is 



W 



-62] + 



uu 



(264) 



(265) 



which since = and $1 = <I>2 reduces to Eq. ([234|). In the general N-body case, the individual wave 
equations become 



i-'l + E ^^^f^]/(E ^) + = (a = 1, N). 



(266) 



For compatibility of these wave equations, we must have 



N 



$a = 



^ab{Xab±), 



where again, P is the total momenturq^ 



N 



a=l 



(267) 



(268) 



The $ in this equation is chosen to have the same functional dependence on S and A as in Eq. ()179p . 
Sazdjian does not deal directly with the potential (hence his free case here with potential added at the 
end) , so part of this work has been to use the potential from the Two-Body Dirac equations in the Sazdjian 



^Note that this dependence of the potential on the part of the potential that depends on component of the relative coordinates 
perpendicular to the momentum of the total system, is allowed as long as the issue of cluster decomposition need not be 
addressed. In that case, where one may separate out a part of the system from the remaining part, it is not meaningful to 
require the potentials to depend on the total momemtum of the total system instead of that of its subconstituents. 
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three-body equations. 



3.2.1 Our Adaptation of Sazdjian's Three-Body Generalization 

In order to apply the work done by Sazdjian to our problem, we have to specialize his A^-body equations 
to the three-body system and derive the appropriate Hamiltonian, eventually ending up with an equation 
that looks very much like a nonrelativistic Schrodinger equation and is thus much easier with which to work. 
He obtains the following approximation for the three-body eigenvalue equation, which is the expanded three 
body version of Eq. (j264l) . Specializing this equation to the three-body case for the baryons and expanding 
terms, we get 

= [e?-™?-3 + 



i(l/ei-Hl/e2 + l/e3) 
ei-mi- 3- 



-2 2 o P2±+^2 



e2(l/£i + l/e2 + l/e3) 

.2 



+4 -ml- 3 f^^ \ A H^i2±,X23±,x,i^), (269) 
£3(^/£l + i-/£2 + i-/£3) 



in which the epsilons, representing the total energy of each quark, are 



W 1 mi — TO2 

^ ^ 3" ^ 3 1 + (w - M)/(mim2(l/mi + l/m2 + l/mg) 

1 TOl — 7713 

3 1 + (ly — M)/(mim3(l/?77i + 1/7712 + ^/ttt-s) 

W 1 7772 ^ "T-l 

^ ~ 3" ^ 3 1 + (77; - M)/(777im2 (1/7771 + 1/^2 + l/mg) 
1 7772 — TT^3 

31 + {w ~ M) / {m2mz{l / mi + l/m2 + I/7773) 

^ , 1 "73 - rni 

^ 3 3 1 + (77;- M)/(mir773 (1/777,1 -f 1/7772-1-1/7773) 

1 7773 — 7772 

3 1 + (tw — Af)/(m2m3(l/mi + l/r)T-2 + I/7773) 
and the potentials $i are linear combinations of the two-body interacting potentials 

^•i = $12(a;i2J_,ei,£2) + $23(a;23_L,£2,e3) 

$2 = $23(a;23±,e2,£3) + *3i(a;3i_L,e3,ei) 

$3 = $3i(a;3iJ-,e3,£i) + *i2(a;i2_L,ei,e2)- (271) 
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This equation is essentially the three-body version of the two-body equation: 



^- — 2^, — + — — ' (^^^^ 



as long as one restricts oneself to confining interactions. 

In the above equations for the epsilons ( Eq. ()270p ). we define 

Pi_L = Pi+Pi- PP 

p 



p = 



P = P1+P2+P3 

E — W — M = Sl + 62 + £3 — TOl — 7712 — W3 (273) 



in order to make them into a more usable form. Through some extensive algebraic manipulations, we obtain 
the much simplified equations of 

El -mi = Efi{E,mi,m2,m3) 
£2-m2 = Ef2{E,mi,m2,m3) 

£3-TO3 = Ef3{E, mi, 1712,1-713) (274) 
in which the //s are collections of functions, given by [5T| 
fi{E,mi,m2,m3) 

Sm^m^ + Sm^m^mi + 'im2m'^mi + mim^E + 2m1m3E + 2m|i?m2 + m^Emi + m3E'^m2 
3 (rniTO2 + 77127713 + 7773777,1 + m^E) {mim2 + 77727773 -|- 77737771 -|- m2E) 



/2(£^, 7771, 7(72, 7773) 

3777^7?727773 + im2m\mi -f m\Em2 + 2m\Emi + im^m^ + 2m3m\E + m^E^mi + m\m2E 

3 (mim2 + TO27)73 -I- 7)737711 + m^E) (7(717712 + 77727)73 + 77737771 + miE) 



/a (-B, 7771,7712,7773) 

3777277737771 -f 3777^77i2'773 -|- 2mimlE -\- m^m^E -\- Sm^ml + m^m^E + 2m\m2E + m2E'^mi 

3 (77I177I2 + 7T127773 -|- m^mi + m2E) (77717772 -|- 77727773 -|- m^mi -\- miE) 



,(275) 
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so that 

ei{E,mi,m2,mz) = mi + Efi{E,mi,m2,m3), 
e2iE, mi, 1712,1713) = 7712 + Ef2iE,mi,7n2,m3), 

ez{E,mi,m2,mz) = ^3 + Ef3{E,mi,7n2,m3). (276) 

Now we want to go back to the full three-body equation fEa. (|269p ) in order to create a Schrodinger-like 
form. Our bound state equation is then rewritten as 

[{el - ml) + 1/82 + l/es) - 3^' + ^1 + *i 



+ (4 - ml) + 1/62 + l/es) - 3 



£1 



£2 

+ {el - ml) (l/£i + l/£2 + l/£3) - 3^^i^^i^]V'(xi2±,X23±,X31±) 

£3 

= 0. (277) 

And once again, in order to simplify down to a Schrodinger-like form, we collect functions into the definitions 

6i(i?, mi,TO2,m3) = (ei(£',TOi,m2,TO3) 

x(l/ei(i;, mi, 7712, 7713) + l/e2(-B, mi,m2,m3) -f l/e3(£', mi, m2, m3))/6, 

b2{E,mi,m2,m3) = {e2{E,mi,m2,m3) + m2) 

x{l/ei{E,mi,m2,m3) + l/e2{E,mi,m2,mz) + l/£3(£', mi, m2, m3))/6, 

b3{E,mi,m2,m3) = (£3(5', mi, m2, m3) + m3) 

x(l/ei(i;, mi,m2,m3) + l/e2{E,mi,m2,m3) + l/£3(£', mi, m2, m3))/6. (278) 



57 



So, our eigenvalue equation is now 



{E[fi{E, mi, 7712, m3)bi{E, mi, m2, ma) 
+f2{E, mi, 7712, m3)b2{E, mi, 1712,1713) 
+f'i{E,mi,7n2,mz)b'i{E,mi,m2,mz)] 

Pl_L + ^12 + ^13 

2ei(i;,mi,m2,m3) 

P2J. + $23 + $ 12 

2e2(-B, mi,m2,m3) 

+ *31 +^"23 11 X 

\}v{xl2^,X2-i^,X3il_) 



2ez{E,mi,m2,mz) 
0, (279) 



and finally, defining 



fi{E, mi,m2, m3)bi{E, mi, m2, ma) 
+f2{E, mi,m2, ms,)b2{E, mi,m2, m^) 

+f3{E,mi,m2,m3)b3{E,mi,m2,m3) = F (280) 



we acquire a (relatively) simple Schrodinger-like form of the Hamiltonian: 

1 / Pi_L + ^12 + ^13 ^ P2± + ^23 + $12 ^ P3_l + ^31 + ^23 , ^281) 

F \2£i(iJ,mi,m2,m3) 2e2(-E, mi, TO2, ma) 2£3(i?, mi, m2, ma) 
As a check, note that for = 



^(0) = /i(0,mi,m2,ma)&i(0, mi,m2,ma) 
+/2(0, mi, m2, m3)62(0, mi, m2, m3) 
+/3(0, mi, m2, m3)6a(0, mi, m2, ma) 

1 (282) 



so that when E ~ 0, this becomes an ordinary non-relativistic Hamiltonian like Ea. (|218p . So now, we have 
gone from Sazdjian's N-body formalism to a condensed three-body one that is much easier to work with 
in the constraint dynamics approach, as the Hamiltonian is now in a familiar form. This equation has 
the distinct advantage that it is like, for the purposes of solving it anyway, a non-relativistic Schrodinger 
equation. It is, of course, still relativistic, but it is now in a form that is much more easily recognizable and 
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usable than the one which Sazdjian left us, Ea. (|266p . It is important to note the recursive nature of this 
equation as this becomes highly relevant in the numerical studies. The $'s are dependent on the e' s and 
the w' s and so we must begin with an initial guess and solve the equation iteratively until an acceptable 
level of convergence is met. 

4 The relativistic three body problem 

It should be noted that this model and our relativistic generalizations below are often referred to as 
the "naive quark model" in that it does not account for the swarm of gluons and quark-anti-quark pair 
interactions directly. Rather, the model has all of the interacting forces existing between each quark-quark 
pair. 

The process of going from a two-body system to a three-body one is not as straightforward as one might 
expect. Since the interacting potentials are limited to each quark-quark pair, there are now three times 
as many terms. Additionally, there are three sets of coordinates (ri — r2,r2 — r3,ri — v^) instead of just 
one distance between the quarks as in the two-body case. This is best treated with a relative coordinate 
substitution to reduce the number of coordinates from three to two. The coordinate transform used in this 
work is similar to and uses the same notation as Isgur's ( 6 ) work, but is not an identical transformation 
due to treating the more general case of all three quarks as different masses, rather than just two. 

The sections that follow describe the methods used in going from a two-body system to a three-body 
one, mostly dealing with the potentials and coordinate transforms. We also describe our Gaussian basis 
functions and the reasoning behind them. Since we are still merely using two quark interaction potentials 
(see Eq. (|183p ). the scalar (5), vector {A) and Darwin ($_d) terms, all simply expand from one term to three, 
to account for all three interactions. However, the spin-spin ($55), spin-orbit {'^sOt^soTt'^sox) and 
tensor {^t) terms require manipulations of total J, total L and total S, that the two-body system neither 
accounts for nor has to deal with and thus have been extensively reworked for the three-body system. 

4.1 Coordinate system transforms and tensor substitutions 

One of the more difficult parts of the three-body problem is also possibly the most fundamental: a 
proper coordinate system. Therefore, this section provides a description of how the coordinate system is set 
up for the relativistic three-body problem. Defining a coordinate system that is both correct and useful for 
a three-body problem is non-trivial and requires extensive manipulations and even some of what would be 
considered mathematical tricks. One of the simplest and most common ways to begin handling a three-body 
system is to redefine the coordinate system so that there are only two relative coordinates instead of three. 
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In the following section we describe the way the coordinate system is defined for the three-body system and 
then reduced to two relative coordinates, plus a "center of mass". We then address how those are further 
simplified with additional coordinate transforms in order to analytically solve as much of the problem as 
possible before going to numerical methods. 

Let us return to Eq. (|28ip . Our goal is to create a coordinate system in which the kinetic terms are 
analytic and the variational principle will be used to solve for the energy eigenvalues. The general form of 
each is 

^ij = 2m^.yS'ij + Sfj + 2ewijAij - A^- + <^di] + L • {cr i+cr j)^ soij + cr^-rijCrj-rijL ■ {cr 1+0-2)^ soTij 
+ a-i-a-j^ssi] + {ia-i-VijCTj -Yij-cri-cT j)(^Ti] + L • (cTj-c^jO^soDy ' (283) 

where 



Sij — 



27 ' 



Stt 6162 
"^'^ ^ "27r.,ln(Xe2+(4^ + i^' ^^^^^ 

and the various $ terms are all functions of Sij^ Aij and their derivatives (explicit forms given in Appendix 
A). Note here how they still account for the asymptotic freedom and linear confinement mentioned earlier. 
The scalar term goes to infinity as r goes to an infinite value, providing confinement, while the logarithm 
in the vector term becomes large at short distance, giving asymptotic freedom (this causes the vector term 
to behave more like ~a/r). We now want to define a coordinate system such that the three coordinates 
Tij become two relative coordinates that can be written in terms of the original Tij distances between each 
quark pair. In direct analogy to what was done in the nonrelativistic case fEq. (|220p . the coordinate system 
is defined as 

p = rs-rs, 

we2 , wez 

(e2+£3)ei (e2+£3)ei 

£3 , 
ri-r2 = ■ p-A, 

£2 +£3 

ri-rg = ■ p + A, 

£2 +£3 



r2 - rs = p. 



£l(£2+£3) ^ £2£3 /„nrN 

, £a = ; , (,2«5j 

w £2 + £3 
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where again, w is the total baryon energy eigenvalue and the epsilons are the individual cm. energies of 
each quark, such that 

w = £i+e2+£, (286) 

and the Cp, and sx can be regarded as reduced energy terms (similar to reduced mass, but energy instead 
since we are writing most of our Hamiltonian in terms of energy). We choose this form (rather than a 
more standard one that Isgur uses) in order to maintain the proper form of our momentum terms from the 
two-body formalism. In this new system, the original Hamiltonian of Eq. (|28ip now becomes 



F 2ep{E,mi,m2,m'i) 2exiE, mi, 7712,1113) 

$12 + ^'13 , $23 + *12 , $31 + *23 



-), (287) 



2ei{E, mi, 1712,1713) 2e2(£', mi, 7712, ma) 2e3{E,mi,m2,m3 

so that the kinetic terms are now simple and easily worked out analytically. Note that this equation is very 
similar to the non-relativistic Hamiltonian Eq. ()218|) in the preceding section. 
We use the variational principle 

i^nim^^m) = {^n\E\^rr.), (288) 

to find the eigenvalues of our Hamiltonian. For greater accuracy, we expand our Hamiltonian basis into the 
matrix equation 

ro=AB^', (289) 

where H and B are both matrices given by 

Bn,n = (290) 



We denote the wavefunction of total J,L and S, where total L is composed of the angular momenta 
associated with the p or A coordinate and total S is composed of the individual spins of each of the three 
quarks, with Si, S2 and 5*3 corresponding to the spins of quarks 1, 2 and 3, respectively, as 

\JL{lph)S{SiS2S3)Mnpnx) = |*„). (291) 
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Therefore a matrix element of the Hamiltonian would be 



{^n\H\-^m) = {JL{lph)S{SiS2Ss)Mnpnx\H\JL'{l'/^)S\S[S'^S's)M'ny^) (292) 

One major advantage of this particular choice of coordinates is that the kinetic terms can be analytically 
evaluated. The matrix elements for the kinetic term 

i_( ^ I ) (293) 

F 2ep{E,mi,m2,m3) 2ex{E,mi,m2,m3) 

are then 



(*„|T|*„) 

sJ. "P^^'"'^'^^ J(n nM(2ip+3)/2('j,,„' ■)(2;A+3)/2r; (i , i ^ r[(2/p + l)/2] 



Fepr[{2lp + 3)/2]y' ' P' ' ^ ^' ' ^2(K + n'p]a2/2)(2ip+i)/2 



2 r[(2^p + 3)/2] r[{2lp + 5)/2] 

V^-^' ^ P2([np + nya2/2)(2ip+3)/2 +"p"P2([np + nya2/2)(2i.+5)/2J 



F£A r[(2/A + 3)/2] V ^ " ^ - ^2([n, + nga2/2)(M)/2 

-r97 |0Wn.2 r[(2^;, + 3)/2] o , r[(2;;, + 5)/2] 

^ ^ + "^^"AaA2([„^ ^ n'Ja2/2)(2/.+3)/2 + "a "a ^^^^^ ^ „'Ja2/2)(2i.+5)/2J> 

X i^'p^'xK"^'>^\^^)^l'A^K"''x Yl i¥x^P^x\LM)Klx^rn,m^- (294) 

The matrix elements for the potentials are not analytic though and must be evaluated numerically. 
However, with the current substitution, the r23 term is relatively simple since r23 = p. The matrix element 
for the spin-independent $23(t'23) potential then is 

(*„(p,A)|$23(p)|*m(p,A)) 

$23 (p) 



(npa2 ) (2;p+3)/2 (^^.a^ ) (2i.+3)/2 / (n'^aS) (2i;,+3)/2 (^/^^2 ^(2i^+3)/2 

r[(2Zp + 3)/2]r[(2/A + 3)/2] V r[(2;;; + 3)/2]r[(2Z';, + 3)/2] 

Y {l,lxmpmx\LM)Y,^''Yl^^ ^ {l%rn'y,\L'M')Y,'r'^Y,fd'pd^X, (295) 



which, due to there being no dependence on the A variable in the potential $23 (p) and performing the angular 
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integrations, reduces down to 



(*„(p,A)|$23(/')|*,„(/'.A)) 
$23 (p) 



(npa2)(2;p+3)/2(„^C,2)(2i,+3)/2 / (n'^a2 ) (2i;,+3)/2 (^,^^2 ) (2;', +3)/2 

r[(2Zp + 3)/2]r[(2ZA + 3)/2] V r[(2z;, + 3)/2]r[(2r;, + 3)/2] 

r[(/A + l'x + 3)/2] (;^+;^+2) -(n„+n:)a^„pV2 

2[(nA + n'Ja2](iA+i^+3)/2^ 



Thus we are left with a function of one variable that can easily be numerically integrated regardless of what 
$23 (p) happens to be. However, this is not the case for the other two terms. The matrix elements of the 
ri2 and ris interactions are, respectively 



(*„(p,A)|$i2(p,A)|*„(p,A)) 

$12(P, A) 



(npa2)(2'p+3)/2(„^a2)(2;A+3)/2 / (n;,a2)(2!;,+3)/2(^^^2)(2;',+3)/2 

r[(2Zp + 3)/2]r[(2;A + 3)/2] V rp;, + 3)/2]r[(2Z^ + 3)/2] 

xp('f+'p)A('^+'A)e-("''+"p)«>'/2-("A+0"jAV2 

^ (Zp/,mpm,|LM)F™''y™^ ^ {l'/xm'y,\LM)Y^'r'''Y;r'^d'pd'\, (297) 



(vI>„(p,A)|$i3(/9,A)|*^(p,A)) 

_ r /(npa2)(2/p+3)/2(„^„2)(2;,+3)/2 /(n'^a2)(2/;,+3)/2(^^^2)(2i',+3)/2 

y$i3(p,A)y r[(2Zp + 3)/2]r[(2/A + 3)/2] V r[(2Z; + 3)/2]r[(2Z^ + 3)/2] 
E (;p^AmpmA|LM)y;;^y,7M3^d3^ ^ {l'/xm'yx\LM)Y,T'^Y^fd^pd'X, (298) 

and so as the potentials are now in terms of two variables, it is much more difficult and time-consuming to 
numerically evaluate this integral. We therefore wish to make another variable change in the ri2 and ri3 
systems in order to rewrite them in terms of a single variable as well. 

What now follows is a brief description of the variable change simplification; more explicit details can be 
found in Appendix D. The variable change used is based on properties of the spherical harmonics and how 
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they relate to spherical tensors. Since spherical harmonics are trigonometric functions, we can get 



pcos6 = z, 

(y,V)'5fp = J-^', (299) 

4 TT 



1/3 1/3 
Ff^ = -W — sin6le""^ = -W — sin6'(cos(^-isin^), 



r cos y sm ( 



(n-vrn-v = ^v^(-+^^)^v^(--y) = J|^(-^+y^)- (300) 

and similarly for the other spherical harmonics. Since we will be using no more than / = 2, we now specialize 
our discussion to ? = 1 with a brief description of / = 2. Therefore, since part of our wavefunction is a 
spherical harmonic (which is a trigonometric function) and a coordinate, the wavefunction can be rewritten 
in spherical tensor form as 

^ M 
= -^n^/4a5/2„5/4^5/2^_„,„J,V2-nxa^V2 ^ (1 1^^,„, |00)p„/„„ (301) 



M 



where 



Pm, = pYr^ip) (302) 

^mx = Ay/"^(A). 

Additional manipulations are still needed in order to work out the expectation values explicitly. For the ri2 
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integration, a new set of variables is defined as 



p' = ri2 = ri - r2, (303) 

, wei we2 
A = —7 ■ rri + — ■ -r2, 

then are rewritten in terms of new primed variables and as tensors, using the same tensor substitution done 
above 

El + £2 

Pm, = —^-P'm,+>^m,, (304) 



, WS2 I £3 

A = -p A , 



(£2 + £3)(£l + £2) £2 + £3 

W£2 , _ £3 

(£2 +£3)(£l +£2)^"^ ^ £2 +£3' 



\ - "'^2 , t3 w 



Note that this is not a new coordinate system but rather a change of integration variables. This means that, 
while we are currently working out the new integral for the ri2 system, we can use a similar substitution 
for the ri3 system and acquire a nearly identical equation with only a few constants changed (constants in 
terms of the integration variable, not overall constants for the full calculation). Before the new substitution 
of the primed coordinates, the expectation value of the potential <&i2(p' ~ T12) is 

(*„|$12(ri2)|«'r«) = (*„|$12(p')l*m) (306) 



16 5 5 / 5/4 5/4 5/4 5/4^ 

^apaA("pi '^Ai "a2 ) 



X 



j ^ (llTOpiTOAi|00)p,*„^^A;,^^ ^ (llmp2mA2|00)p„^^A„^,, 



mpimxi mp2rnx2 
Xe-Ki+'V2)">'-("Al+nA2)aiA^^3^^3;^^^^(^/)^ 

And so now, using the derived relationships between the primed and un-primed coordinates p and AfEa l302l303p . 
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the expectation value becomes 



(*n|$12(p')l*m) 



where for simpUcity 



-'^ — g^'^p'^xK^pl "■Al "p2 )' 



ei+£2 2 (£2 +e3)(ei +£2) 

(n^i + np2)ap , £1 ^ , . , 2 ^8283 
7: ( \ ) + {n\i + nx2)a 



ei+e2 (£2 + e3)^(ei +£2)' 

__^3_-^ 
£2 + £3 ' 



^ ^ (npi+np2)ap ^ {nxi + nx2)al ^ £3 ^2_ (309) 



However, there is still a troublesome A' • p' cross term in the Gaussian wavefunction itself. In order to 
eliminate this, one final coordinate change is used 

V = x-ip', 

= Xm-^p'm- (310) 

The explicit details can be found in Appendix D; the end result is 



66 




+ 2^, 



) 



]}) 



J 4 16 

X (llmpimAi|00)[ ((ll"ip2mA2|00) 




X 



Y -(llmA2mp2|M(1100|/0) 



( 







1 1 I 



) 



( 



mpi mAi m 



1 



1 



) 



)] (311) 



So now the expectation value is in terms of just one variable, which can be numerically evaluated easily. 
The ri3 term simply requires changing the a, b, and c constants and the epsilons and the term is trivial 
since the coordinate system is defined in terms of r23. Again, these are listed in Appendix A. The potential is 
now in terms of just one variable, so almost regardless of what potential is used, the numerical calculations 
will be fairly straightforward. Thus, the coordinate system has been defined and transformed in such a 
way as to make a good deal of the problem analytic while keeping what is not analytic still easy to solve 
numerically. 

We can follow a similar method for Z = 2, by constructing new states as functions of the I = 1 states. 
The explicit forms of the spherical harmonics for Z = and I = 1 are 




(312) 



which can be used to rewrite the / = 2 spherical harmonics as 
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Y- 



Y: 



^2 



^2 



Y- 



■ sm f cos fe 



i^f (Scos^e- 1) = (yO)2-Ly^_ (^^0)2^5^^ 



1 

2 V 27r 
1 ^15 



(313) 



So then, since our wavefunction for I p — l\ — 2 \& of the form 



I^«(P,A)) = 



ny4«y2„V4^7/2^2^2g-„,a>V2-n.aiAV2 



^ (22mpmA | LM) Is"" ^2"^ > 



ITlplTlA 



we can use the original tensor substitution of Eq. p02p to define 



(314) 



P -'2 
P -'2 



\Y,rv2jf^v2^r-;^pi 



(315) 



Therefore, we can follow the same procedure outlined for 1 = 1, the only major difference being that the 
algebra is more complicated. However, since we're going to have to do at least one numerical integration 
anyway, we can simply use a program such as Mathematica (or Wolframalpha) to analytically integrate 
what is analytically integrable and then numerically integrate the rest without doing an absurd amount of 
algebra. This method should actually hold to allow for higher order I states to be done relatively simply, 
since essentially all that will become more complicated is the algebra arising from the substitutions (this 
may be part of future work on this subject). 

With the matrix elements defined for a general potential and for analytic kinetic terms, we now need to 
explicitly define our potential model. 
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5 Three Body Potential 

Conceptually speaking, the approach one would take to go from a two-body system with the formalism 
we have described to a three body one is straightforward. The problem is now treated as three two-body 
problems, so the overall form of the potentials remain mostly the same. However, the explicit calculations 
become considerably more complicated since there are three different sets of couplings to consider, three 
different sets of coordinates (which were reduced to two in the previous section), and three different operators 
for each potential term. This not only increases the magnitude of the calculations involved, it requires 
entirely new approaches as some methods that work very well for the two-body system can no longer be 
used. As everything is now in terms of the r23 variable, the other terms must be appropriately recoupled 
to properly use the operators. The two-body potential derived by Crater [5] can be written as 

$ = 2m^S + S'^ + 2e.^A - + <i>D + ■ (cri+cr2)$so + cri-fcra-fL ■ (cri+cr2)$sOT 

-I- cri cr2^ss + (3cri-fcr2 ■r-(Ti (T2)^T + L • {(Ti-cr2)^gQjj + il,-(TiX(T2^sox , (316) 

where the ^ terms are all functions of the vector potential, scalar potential, and quark masses (these are all 
defined explicitly in Appendix A). The three-body potential is of similar form and essentially just triples 
the number of terms, so instead of having 1 and 2 there are interactions for all three quarks. For the 
relatively simple vector (A), scalar (S) and Darwin (<&£>) terms this is almost trivial, as there are no direct 
spin-dependent operators; however, the spin-spin, tensor, and spin-orbit terms require extensive reworking, 
which are described in the following sections. 

5.1 State Couplings and Operator Methods 

We will now describe how we set up our three-body states in such a way as to use the potential operators. 
In order to simplify our numerical calculations, it is helpful to note that the potential terms are products 
of a term involving the coupled angular momentum operators and a coordinate dependent terms that have 
trivial operator dependence . This allows us to use the operator angular momentum operators on a specified 
state and just get a number back that depends on the angular components of the state itself and not any 
radial components, so that the numerical integral itself does not involve any angular momentum operators. 
Thus, our potential terms separated into operator and non-operator pieces are given in table 2 and and the 
explicit forms of the <i> terms are given in Appendix A. For the purposes of working out the operators, it is 
sufficient to simply consider them to be radial functions that are unaffected by the operators. 
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Table 2: Potential terms, operators and non-operator components 



Potential Term 


Operator component 


Non-operator component 


Spin-Spin 


(TiCT j 


3'ss(r.y) 


Spin-Orbit 


Lij-{(Ti+cr j) 




Spin-Orbit Difference 


L,ij-{(Ti — cr j) 




Tensor 


icn-YijCTj -Tij— (Ti-cr j 




Spin-Orbit Cross 


%\jij - (T i^ (T j 




Spin-Orbit Tensor 


o'i'^ijcrj'^ij'i-'ij'icr i+cr j) 


$soT(r.y) 



Table 3: Spin-flavor wavefunctions for the baryon spectrum 



N 


J 


L 


S 


SF 




8 


i 

2 





i 

2 


^(fx' + '/'"x")V'o 




10 


3 
2 





3 
2 




*2 


8 


i 3 


1 


i 


m'x"+ci^"x')f + {<!>' x'-rx"W] 


*3(J=i),*4(J=f) 


8 


2 ' 2 ' 2 


1 


2 


-l^lcjj'x'ip' + (j)"x'tp"] 


*5(J-|),*6(J=i),*7(J-|) 


10 


1 3 

2 ' 2 


1 


i 

2 




*8(J=i),*9(J=§) 


1 


1 3 

2 ' 2 


1 


1 
2 




1'io(J= i),*ii(J= 1) 



For the baryons we have considered, there are a total of eleven different wavefunctions, which represent 
all possible spin-flavor couplings for the various particles. The form of these is given in table 3 where it 
has been split into three components: spin, flavor and space (represented by x, 0, ip, respectively), explicitly 
defined in Appendix B. As there are six interactions to consider and three couplings per interaction (we are 
using two-body operators, so there is a 1-2, 1-3, and 2-3 term for each operator), there are a total of 198 
possible interactions to consider. Fortunately, many of these are similar or trivial and so the number that 
must actually be worked out explicitly drops considerably, but there still are a quite a large number that are 
non-trivial. The eleven wavefunctions are given by [4 7) where ipQ is a total L — Ip = l\ = wavefunction 
and ip' and ip" are total L = 1 and Ip = 1 oi l\ = 1 wavefunctions, respectively {L = 1 state have parity of 
— 1, so Ip and l\ cannot both be 1). In addition, 0',(/>", 0^ and 0" are all purely flavor wavefunctions and 
X** s-re all purely spin wavefunctions, having total S ~ 1/2, 1/2 and 3/2, respectively. These merely 
contain all possible combinations of flavor or spin so that the product of the two gives all possible spin-flavor 
couplings. All of these wavefunctions are orthogonal to the others in the set (that is, x' is orthogonal to 
x" and x'' etc.). Explicit forms of these terms are given in Appendix B. The quark flavor combination of 
the (p terms is different for each baryon, but since the operators we use do not affect the flavor, it does not 
matter what they are explicitly for the purposes of this calculation. 

There are two methods we use to determine the affect of these operators. One is a simple ladder operator 
approach and the other involves use of the Wigner 9j recoupling coefficients. Both methods are always valid, 
but not necessarily always useful due to how the operator form affects each individual wavefunction for the 
ladder operators. It is worth noting that having two methods be viable also allows for a good check. The 
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ladder operator method works out simply for everything but the Spin-Orbit Cross term (Eq. ((TH7])) due to 
how we can have a wavefunction independent of total M . Thus, we can set M — J and force Mg = S. For 
the states ^E"!, ^'2, 5'4^'7, ^"9 and ^fn, this means that any operator that changes total Ms will be orthogonal 
to the original wavefunction and thus we can eliminate any term that does change total Ms- As most 
readers are familiar with ladder operators, we relegate most of the details of this method to Appendix C. 
The other method (9j coupling) though is not as familiar and so we will describe it in detail. Due to having 
matrix elements for the two-body problem already defined by Crater [?] , the difficult part of this problem 
is recoupling the state into one which can use these matrix elements, since it is not readily clear how they 
operate on a three-body problem. 

5.2 Spin-Spin 

This section deals with how the spin-spin couplings are defined and worked out. As we will see in the 
spin-orbit term (and others as well), the use of 6j recoupling coefficients is necessary. The default state is a 
2-3 coupling one due to how we have defined our coordinate system, so in order to work out the spin states 
for the other couplings (1-2 and 1-3) the state must be changed to reflect those couplings, since while the 
S2 • S3 term simply yields a number, the Si • S2 and Si • S3 do not. The recoupling is done using 6j (Racah) 
coefRcients. Here, we work out an example state (the proton) to show how this is done. Full calculations 
and results are in Appendix B. 

Note on notation: the general quark spin state is 

[S1{S2S3}^-]^ (317) 

where Si represents each quark (1,2 or 3), 6*23 is the coupled spin of quark 2 and quark 3, and S is the total 
spin of the system. The curly brackets denote which quarks are sub-coupled (that is, which two are coupled 
first to a Sij coupling before being coupled to total S). For example, 

[S1{S2S3}'^1'' (318) 

would indicate that quarks 2 and 3 are coupled to S'23 and quark 1 couples to that state for total S, while 

[{S1S3}^^^S2]^ (319) 

would indicate that quarks 1 and 3 are coupled to Si^ and quark 2 couples to that state for total S. 
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The spin dependence of the wavefunction is contained entirely within the x'lX" and x'^ terms. These 
X spin terms couple with a flavor term we call (j) in a, unique way for each set of baryons that share similar 
properties. Explicit forms of these are given in Appendix B. Therefore, all we are concerned with for the 
spin-spin operator are these x terms. The various spin terms are written as 



X'ilTlS 

x"("^s 
X"{ms 

X^{fns 
X^{ms 



= ^) = ^(nt-;n), 

= -^) = ^(m-m), 

= i) = -L(2tn-4n-nt), 

= -i) = -L(-2ut + m + m), 

= 3/2) =tn, 

= i) = -L(;tt + nt+m), 

= -^) = ^(ut + m + m). 



(320) 



These can be written in a Clebsch-Gordan sum form for the 1-2 coupling as 



X'{fns 
x!'{rns 

X^i'ms 



^) = x'(ms = -i) = [{W°S3]^/^ 
l) = x"{ms = -l) = [{s^S2yss]'/' 

\) = X'{ms = \)= X%ms = = {{s^s^Vs^f^^ 



(321) 



The easiest way now to convert these to 1-3 and 2-3 systems is to use a 6j recoupling coefficient, which is 
defined as 



|[{ac}^-6]^) = ^(-l)^-^+Jc+J..+J.c^(2j^, + i)(2J„, + 1 




\[{aby^''cy) 



(322) 



where 



Ja Jb ^ab 



(323) 



7 



y J jc Ja 

is a 6j symbol, similar to 9j and 3j symbols. Like the 9j symbol, the 6j allows us to recouple a state into 
a desired form, which allows the operators to work simply. The diff'erence between it and the 9j is simply 
that the 6j recouples three different angular momenta while the 9j recouples four (and in case the reader 
is curious, there are 12j symbols that change couplings among 5 angular momenta and 15j that change 
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couplings among 6, but those are not useful to us here). Thus, to go from 1-2 to 1-3, we have 



X = 



[{S^s,rs,]'/^ = ^(-l)l+^-+V(2^13 + l)(l) I ^ ^ [{«1^3}""52]^/^ 



X = 



i i 



1 1 o 

2 2 -^IS 



2 2 



X = 



[{SlS2}^S3]^/' = ^(-I)l+S- + V(25l3 + 1)(3) M J [{S1S3}^"S2]^/' 
^[{.1.3}°S2]^/^-i[{.1.3}^.2]^/^ 

[{S1S2}^.3]^/^ = V(-l)l+^- + V(2^13 + l)(3) I [{51.3}^- .2]^/^ 

Si3 \3/2 1/2 1 



(324) 



and similarly for 1-2 to 2-3, we have 



X' = [{sis,rss]'/' = E(-1)'^'^^^V(2^23 + 1)(1) 



i[W3}%]^/^ + ^[{32.3}^.l]^/^ 



1 1 

2 2 



2 2 



1 1 



'23 



[{8,82}' 8s]'/^ = ^(-l)l+^- + V(2523 + l)(3) r ' [{s2S3}^-Si]l/2 

^23 \ T o 1 / 

%{828,Y8,Y'--]-[{S283Y8,Y'^ 



X 



= [{8l82}'8sr/' = ^(-l)l+^- + V(2523 + l)(3) 

5^23 

= [{8283}'s,f/^ 



^1/2 1/2 523^ 



3/2 1/2 1 



[{s2S3}^-si]'/' (325) 



So now that spins i and j are coupled, the operator Sj • Sj can operate simply and just gives a number 



_ 02 _ c2 n — 3 _ 3 

C2 _ c2 _ c2 9 _ 3 _ 3 



1 



(326) 
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Also note that in terms of ladder operators, 



S, • Sj — —{Si-^-Sj- + Si-Sj-^-) + SizSj, 



•Z-i 



(327) 



which works out simply for every ^„ listed above. The 6j recoupling is then not necessary for the spin-spin 
interaction, but it does become necessary in the others. Details and tables of all spin matrix elements are 
given in Appendix C. The spin-spin is, unfortunately, the simplest of our interactions. We have described 
two possible ways of obtaining the operator matrix elements here and now we continue on to the spin-orbit 
interaction. 

5.3 9j Couplings 

Here will we describe how the 9j recoupling coefficient is used for our three-body problem. We begin with 
a totally general state with four angular momenta of 



where ji and ^4 couple to jV, 3i and couple to jg and finally jV and couple to a total jgmg state. If 
instead we wanted to couple the state to 



then we can use a 9j recoupling coefficient. The definition of the transformation is (essentially a big 
Clebsch-Gordan coefficient) 



I {{jm)j7, {32k)k)39mQ) 



I((ilj2)i3, (j4i5)j6)j9TO9) 



I((jlj4)j7, (J2j5)j8)j9m9) 



X] X] l((.?'iJ2)i3, {kjh)k)h'mc,) {{{3132)33, {k3b)k)hmQ\{{hk)37, {323b)3&)39mQ) , (328) 



h je, 



where, like a CG coefficient 



{{{hh)hAkk)k)krnQ\{{3ik)37,{32k)k)krnQ) 



r 



31 32 33 



V(2i3 + l)(2j6 + l)(2i7 + l)(2i8 + 1) < 



34 35 36 



> . 



(329) 



37 3s 39 
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Note that this 9j coefficient is the reason we have written our original state as |((iii4)j7, (j2i5)j8)j9"^9) 
instead of |((jij2)j3, (j4j5)i6)i9'7i9), it allows greater clarity in the 9j symbol itself due to how the numbers 
read simply left to right (that is, ji, 32, hi ji--)- 

So now, since our initial state is the total JLM state 

\{{S,S,k)S,{klj)L)JM), 

and we want a state of total J coupled in such a way that the operators derived by Crater et al. i4j work 
simply 

\{{S,li).h,{Sjklj)Jjk)JM), 
we can use the 9j symbol to recouple the states as 



\{{S,Sju)S,{klj)L)JM) 

J2 \iiS^k)J^, {S,kli)Jjk)JM){{{S,h).hi {Sjuh)Jjk)JM\{{S,Sjk)S, {l,lj)L)JM) 



Ji Jjk 



X < 



Si 1% Ji 

S L J 



> \{{S,h)J„iSjklj)Jjk)JM) 



(330) 



where < 



Sjk Jjk 

S L J 



> is a 9j symbol. For example, if our state was a 2-3 coupling, then the 9j recoupling 



would look like 



\{{SlS23)S, {hlp)L)JM) = |((5i/a)Ji, iS23lp)J23)JM){{{Silx)Jl, (523^p) J23) JM|((5i523)5, {lph)L)JM) 

Jl .J23 

= E v/(2^i + 1)(2J23 + 1)(2^ + 1)(2L + 1) 



X < 



S Ix Jl 

S23 Ip Jjk 

S L J 



} \{{Sih)Ji,{S23lp)J23)JM). 



(331) 
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Every row and column of the 9j symbol itself must conform to the triangle condition as well (angular 
momentum Sudoku, basically), so while it appears to be a complicated and intimidating method at first 
glance, it quickly reduces down to something very manageable. It is important also to note that one of the 
advantages to using 9j symbols as opposed to combinations of 3j symbols is that the 9j is independent of 
any M values. Explicit calculations and tables of the 9j symbols are given in Appendix C. 

5.4 Spin Orbit 

This section defines and describes the spin orbit coupling for the three-body system. The focus is on 
how to handle the L ■ S coupling for the three-body system, as the other terms are simple variables. This 
is where the 9j symbols begin to come in very handy. In order to use the operator Lij ■ Sij in a simple way, 
we need the wavefunction to be in a total Jij — Lij + Sij state. Then, the spin-orbit operator working on 
a state {Uj Sij Jij) gives 

hj • SijlhjSijJij) — [Jij{Jij + 1) ^ hj{lij + 1) ^ Sij{Sij + l)]\lijSijJij) 
lij ' Sij\lijOJij) — 

hj ' ^ijl^^ij^ij) ^ (332) 



Since we have three couplings, we have three different hj terms. Due to how we are completely separating 
the operators from the spatial components of the wavefunction, these states can still work out relatively 
simply for the 1-3 and 1-2 couplings (recall that our "default" coupling is 2-3, defined in terms of the p 
variable). The momenta associated with these coordinate changes are, by definition, the momenta that are 
conjugate to each of them. So then, our three angular momentum terms would be (using the definitions 
given in Ea l303|) 



Ll2 


= ri2 X pi2 




= f 13 X Pl3 


L23 


= ?'23 X P23 



(333) 



Our substitutions for the 1-2 and 1-3 systems give a linear combination of wavefunctions-as the coordinate 
transform itself is a simple linear combination-of either Z = 1 or ^ = 0, so the calculations for those states for 
the spin-orbit ends up being the same as that for the 2-3, just with a few more terms. 

We can then apply the 9j symbols with the spin-orbit operator to obtain a value for the operator on a 
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given wavefunction. For a state of Ijk = or Sjk = 0, the operator gives zero, which greatly simphfies 
calculations. Therefore, the spin-orbit operator Ijk-Sjk working on a state \{{SiSjk)S, {lilj)L)JM) gives 



\,u-^,k\{{S,Sju)S,{hl,)L)JM) 



^ ^ Ijfc • S,fc ^(2 J, + l)(2J,fe + 1){2S + 1)(2L + 1) 



Ji Jjk 



\ 



Ij Jjk 



} \{{Sdi)JrAS,kl,)Jjk)JM), 



(334) 



5* 



L 



J 



where we have used the 9j recoupling coefBcients to change the state into one which the spin-orbit operator 
is merely a number. Details of these calculations are given in Appendix C. 

As explained earlier, some of the wavefunctions allow the use of ladder operators on the base functions 
without any concern to couplings. Where this is doable it is considerably easier for the spin-orbit state since 
when the 5+ or S- matrix elements are zero, the spin-orbit operator reduces to 



or just the z components of the operators. Again, both methods are always viable, but the underlying 
structure of the wavefunction determines which is more useful for a given calculation. 

The details of the ladder method are also given in appendix C. We have recoupled the state for a general 
wavefunction into one which the spin-orbit operator is simply a number by using the 9j symbols. Now that 
the two simplest terms are completed, we move on to the more complicated tensor interaction. 

5.5 Tensor(<l>T) 

The so-called tensor term, [3cri-fcr2 -f — cri-cr2]$T term involves couplings between different possible 
L states. Since our potentials are still two-body interactions, the effect of the Srij = Scn-^ijCTj ■^ij—o'i-o'j 
operator remains the same for each interacting pair. Like the spin-orbit term, it involves total S and total 
L of the two-body system, so that in order to use the matrix elements from a two-body system derived by 
Crater et al. ([TJ [21 [3]), we need to recouple it using the 9j symbols as well. From ([H]), we have the 
following definitions for the Sxij operator working on an llsjm) state, where l,s^j, and m are the two-body 
orbital angular momentum, spin angular momentum, total angular momentum and z component of the total 
angular momentum, respectively: 



Ijk ■ Sjk — Ijk ■ {Sj + Sk) — ljkz{Sjz + Skz) 



(335) 
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ST\lOjm) = 0, 



Srljljin) = 2\jjlm), 



(336) 
(337) 



StIj + Ujm) 



2j + l' 



2j + l 



ST\j-njm) = -^_i|j-lljm) + |i + llim)6^^|^. 



Applying the tensor operator after using the 9j recoupUng coefficients then gives 



(338) 



[3<jj-r(Tk t - (Tj-<JkMSiSjk)S, {lilj)L)JM) 



^^[3<7,-f<Tfc -f - aj-ak]yJ{2Ji + l)(2J,fe + 1)(25 + l)(2i + 1) < 



x\{{Sih)Ji,{Sjklj)Jjk)JM). 



Si li Ji 

S L J 



(339) 



Explicit forms of this equation change depending on which baryon state to which it is applied, so the details 
are again in Appendix C. 

The ladder operator is also used for some states where it is useful. The full form of the operator is 



3Si-rSj-r 



3( 



Si-y. + Si- Sj-y. + S. 



3- ^2 



Si-\- + Si— Sj-\. Sj 



Si+ + Si. 



Si+ — Si- Sj-y + Sj- 

+ o + 



2 2i " 2 

Si-y — Si- Sj-y — Sj- 2 I Si-y — Si- 



+ 



2i 2 

^ >~>izl z' X T "Jj; 



V ' 



Sj+ — Sj- 

2i 



2i 2i y 2i 

'^zfy + SizTzSjzTz) — Sj-Sj, 



'Sjz'^x'^z 

-SjzTyVz 



(340) 



which when we only consider terms that do not change total Ms becomes 



Si+Sj- + Si-Sj+ ^.^2 ^ ^ ^^^^^^^ ^^^2 Q _ S..S . 



(341) 



that, as one may expect from, becomes simple in most cases where it can be applied because sin^ 6' + cos^ 9 = 
l(and thus where this is useful the 6 does not matter since it disappears from the equation). This dramatic 
simplification is only possible in the highest possible J states for each set of baryons, as when J is at a 
maximum value, we can set M = J (since everything is independent of total M), which in turn forces total 
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Ms to be its highest possible value, meaning that any change to total Ms on the ket would cause it to be 
orthogonal to the bra and therefore those terms would not contribute. 

Thus, the adaptation methods of the two-body tensor operator to a three-body system is now well-defined. 
We continue on to the fourth and fifth potential operators, the spin-orbit difference and the spin-orbit cross. 



5.6 Spin-Orbit Difference and Spin-Orbit Cross {^sod.^sox) 

In a similar manner to the tensor term, here we describe how the spin-orbit difference and spin-orbit 
cross operators work on the total {Isjm] state. From Eq. (|186[) we have that 



L- (rTi-rTj)$50D|jm) = \jlln)[{l-{j-lf){{j + l){j+l + 2))]^sOD 
L • [cTi-a^)^soD\]lln) = \jlQn)[{l - {j - + + I + 2)}]^sod 



(342) 



tL-a,xa,<^SOx\jlOn) = |j7l7i)[{l - (j - lf}{{j + + I + 2)}]$sox, 



ih-(TiXcrj^sox\jlln) 



-\jlOn)[{l - [j - 1)'}{U + l)U + l + 2)}]$ 



SOX- 



(343) 



As can be seen from the operator equations, these states will account for spin mixing effects, coupling spin 
zero states to spin one states. Recoupling using the 9j symbol gives 



L • (<Tj-<Tk)|((5,^,fc)5, iklj)L)JM) 



^ • {aj-cT^)^{2J, + l){2J,k + l)i2S + 1)(2L + 1) < 

Ji Jjk 

x\{{S,k)J,,{Sjklj)Jjk)JM). 



s L J 



(344) 



These particular states also lend themselves very well to the ladder operator method, much like the spin-orbit 
state. For the spin orbit difference, the only non-zero operator (that doesn't change total Ms) is 



L • (cTj — CTk) =Lz{a jz — CTfcz), 



(345) 



and for the spin-orbit cross term we also get 



(346) 
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These ladder operators are however only usable on about half of the possible states. Details of these and 
the other operators are given in Appendix C. 

The spin-orbit cross and spin-orbit difference terms have been described from the purely two-body for- 
malism to the three-body formalism, and so now we continue to the last term, the spin-orbit tensor term. 

5.7 Spin-Orbit Tensor 

The final piece of our potential is the spin-orbit tensor term,cri-f cr2-fL • (cri+cr2). Unlike the other 
terms, this one does not lend itself well to the ladder approach, so all of the states have been worked out 
using the 9j method. The operators on an \lsjm) state are ( [49j ) 



orfcTj-fL • {cr-+(Tj)\jljm) 
a-i-rcr j-rL ■ {cr^+crj)\j ~ lljm) 



[j{j + l)-l{l + l)-s{s + l)]\jljm), 
+ 1)- 1(1 + l)-s{s + l)] 



o-.-foTj-fL- (cr^+crj)|j + lljm) = + 1) - + 1) - s(s + 1)] 



2 j + 1 ' 



2j + l 



(347) 



As with all of the potential terms, the individual values of these operators will vary, but we can write a 
completely general state as 



J2 '^j-f'^fc-fL • {cT^+(Tk)^J{2J, + 1)(2 J,fc + 1){2S + 1)(2L + 1) < 

Ji Jjk 

x\{{S,k)J,,iS,klj)Jjk)JM). 



Si 1 1 Ji 
S L J 



(348) 



Details of all calculations for these operators are in Appendix C. 

This completes our derivation of the individual Hamiltonian matrix elements. We have written the 
three-body problem in terms of two relative coordinates and shown how they can be transformed for each 
interacting pair. This allows a description of the methods used to adapt the ([Tl[31[3]) two-body potential 
operators of Crater et al. derived for the meson spectrum to the three-body problem. We have used 9j 
symbols to recouple the wavefunctions into usable form and also have used ladder operators for some of the 
wavefunctions (the ones on which it works more simply than the 9j method). It now falls to numerical 
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calculations to obtain a fit for the spectrum using these equations. 

6 Numerical Results and Comments 

The expectation value of the Hamiltonian cannot be evaluated analytically, so it falls to numerical 
studies to acquire an explicit number. We use a Monte Carlo approach combined with a simple gradient 
method to obtain a best-fit chi-squared for the spectrum as a whole, as compared to current experimental 
data. The following sections describe the numerical methods used and give the results after using said 
methods. 



6.1 Methods and Pcirameter Values 

The numerical calculations were done using a Monte Carlo approach followed by a gradient method to obtain 
a least square fit for the spectrum as a whole. We originally attempted to use a more simplified gradient 
approach but it quickly became apparent that the functions are far too sensitive to changes and thus will 
get "stuck" in a local minimum much too easily without some other approach. So, we adopted a Monte 
Carlo routine that would trigger whenever the gradient approach found a new best fit in order to ensure 
we were reaching the best results for our theory. The integrations were done numerically using Gaussian 
Quadrature and the parameters ap and a\ were minimized by the Nelder-Mead simplex method, though it 
is worth noting that the a parameters do not generally vary much from the analytic result if one were to 
use a harmonic oscillator model and also that as the size of our matrix increases, the actual value of these 
parameters does not affect the fit as much, becoming irrelevant at an infinitely large matrix. Due to how 
we use the additional basis functions to essentially spread out the wavefunction, the results given in tables 
4,5,6 and 7 are for N = 3, which corresponds to a 26 x 26 matrix, as beyond this we observed the overall 
energy and thus the fit would change very little, if at all. 

Our model has a total of 8 parameters, given along with the values for Crater et al. where u,d,s,c and b 
are the masses of the up, down, strange, charm and bottom quarks, respectively. A, K and B are coupling 
constants in our model. It is worth noting that our model has significantly fewer parameters than most 
models, with only 8 total and 5 of those being universal to any model (the quark masses themselves). This 
is significant simply because of the way that a fitting routine will work; the more parameters one has, the 
easier it is to fit the data. Or as the mathematician John von Neumann put it: "with four parameters I can 
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Table 4: Parameter values 





This work 




Reference 2 


u 


166.6 (McV) 


u 


55.7 (MeV) 


d 


165.4 (McV) 


d 


55.3 (MeV) 


s 


337.5 (MeV) 


s 


249.9 (MeV) 


A 


285.8 (MeV) 


A 


421.8 (MeV) 


c 


2050.3 (MeV) 


c 


1.476 (MeV) 


b 


5301.5 (MeV) 


b 


4.844(MeV) 


K 


18.1 


K 


4.198 


B 


100.6 


B 


0.05081 



Table 5: Low lying baryon states 



Baryon 


j 


Theoretical Mass (MeV) 


Experimental Mass(MeV) 


Exp-Theory(MeV) 


P 


1/2 


961 


938 


-23 


n 


1/2 


962 


939 


-22 


S+ 


1/2 


1259 


1189 


-70 




1/2 


1261 


1192 


-69 




1/2 


1262 


1197 


-65 




1/2 


1388 


1314 


-74 




1/2 


1389 


1321 


-68 


A" 


1/2 


1093 


1125 


32 


A++ 


3/2 


1266 


1232 


-34 


A+ 


3/2 


1266 


1232 


-34 


A'' 


3/2 


1267 


1232 


-35 


A- 


3/2 


1268 


1232 


-36 


S+(1390) 


3/2 


1393 


1383 


-10 


E"(1390) 


3/2 


1394 


1384 


-10 


I]"(1390) 


3/2 


1395 


1387 


-8 


2^^(1530) 


3/2 


1505 


1531 


26 


S-(1530) 


3/2 


1506 


1535 


29 


n- 


3/2 


1609 


1672 


63 



10 



The model of course, is still expected to 



fit an elephant, and with five I can make him wiggle his trunk. "[I 
be accurate regardless of the number of parameters, but it is worth noting in this work. In addition, there 
are only two parametric functions that define our potential model, so this further reduces the effect that a 
large number of parameters would have on our model. 



6.2 Results and Comparison to Experiment 

The complete results of our model are given in tables 5-7. As the purpose of this work is to test if the model 
used in the two-body case works well for the three-body, we are only using those baryons which have a three 
or four star rating by the Particle Data Group, meaning that they are fairly well-known. 

The ground-state baryons are generally slightly high energy-wise for the first 8 and this is most likely to 

^"^It should be mentioned though that due to the way our potentials are defined, simply increasing the number of parameters 
should not give a drastic increase in the quality of the fit 
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Table 6: Higher order baryons states 



B&ryoii 


J 


Theoretical Mass (MeV) 


Experimental Mass(MeV) 


Exp-Tlicory 


N(1440) 


1/2 


1508 


1420-1470 


-63 


A(1600) 




1684 


1560-1700 


-54 


Sfl660) 




1660 


1630-1690 





S(1690) 


1/2 


1787 


1680-1700 


-97 


Afieoo) 

i 1 \ X v7 v_/ v_/ / 


1/2 


1535 


1550-1700 


90 


I](1670) 


3/2 


1679 


1665-1685 


-4 


\ J- 'J 'J 1 


1/2 


1561 


1525-1545 


-26 


A(1670) 


1/2 


1678 


1660-1680 


-8 




1/'^ 


1674 


r73()-18()() 


91 


S(1775) 


5/2 


1673 


1770-1780 


102 


Nfl520') 

±1 y _i_ \j \j 1 


3/2 


1562 


1515-1525 


-42 




3/2 


1677 


1685-1695 


13 


S(1820) 


3/2 


1780 


1818-1828 


43 


NCieso) 


1/2 


1549 


1645-1670 


108 


A(1800) 


1/2 


1667 


1720-1850 


118 




1 12 


1679 


1800-1960 


201 


N(1700) 


3/2 


1581 


1650-1750 


119 


N(1675) 


5/2 


1561 


1670-1680 


114 


A(1830) 


5/2 


1690 


1810-1830 


130 


S(1950) 


5/2 


1791 


1935-1965 


159 


A(1620) 


1/2 


1554 


1600-1660 


76 


A(1700) 


3/2 


1549 


1670-1750 


161 


A(1405) 


1/2 


1410 


1402-1410 


-4 


A(1520) 


3/2 


1680 


1518-1521 


-160 



allow the following 10 to be fit relatively accurately. This is not a surprising result of our model due to the 

fact that since wc arc using no purely 3-body potentials, the only difference between these sets of baryons is 
the spin-spin interaction. 

As for the higher order baryons, these are not as well-known and the experimental values are only given 
in a (sometimes wide) range, the difference between the experimental and theoretical values is given as the 
difference between the theoretical value and the average experimental value. This data is given in table 6. 

The higher order baryons fall within an acceptable range on the whole, though there are a few outliers. 
Of important note is that our model does fit very well the often troublesome A(1405) particle. The other A 
particles however are, as before, missing some sort of interaction that will aid in differentiating among them. 

In addition, we fit a the well-known charmed and bottom baryons, given in table 7. These agree relatively 
well with experimental data. 



6.3 General Comments 

The purpose of this work was to see if the two-body model used for the meson spectrum can be purely 
applied to a three-body problem while remaining in what is known as the "naive quark model," that is. 
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Table 7: Charmed and bottom baryons 



Bciryon 


J 


Theoretical Mass (MeV) 


Experimental ]VIass(]VIcV) 


Exp-Theory(MeV) 




1/2 


2397 


2454 


57 




1/2 


2411 


2520 


109 


A^f2286~) 




2402 


2286 


-116 


Ac(2595) 




2428 


2595 


-167 


A^f2880~) 


1 12 


2593 


2880 


-287 


S.f2467) 


1/2 


2400 


2467 


67 


S^f2470') 


1 12 


2401 


2470 


69 




1 12 


2539 


2645 


106 


Sc(2790) 


1/2 


2672 


2790 


118 


Sc(2815) 


1/2 


2678 


2815 


137 


f7c(2695) 


1/2 


2550 


2695 


145 


f7c(2770) 


1/2 


2567 


2770 


203 


S]fc(5829) 


1/2 


5832 


5832 





Sfc(5836) 


1/2 


5833 


5836 


3 


Sfc(5790) 


1/2 


5834 


5790 


-44 


17fc(6071) 


1/2 


5845 


6071 


226 



without considering fully three-body forces, ignoring any effects on the spectrum due to the instability of 
the particles and focusing only on interactions between one pair of quarks at a time. We have shown that 
this is indeed possible as our fit was relatively good compared to other approaches( [9] ) ( ) f [52] ) ( [53j ) and 
to the experimental data itself. Furthermore, the discrepancy between the theoretical and experimental 
results can be explained (qualitatively at least) as one overall missing piece that skews the fitting routine. 
We believe this to be the reason that the quark masses are higher than those done by Crater et al. [39] In 
general, though the initial octet is high, the splittings are correct among them. The masses of the up and 
down quark themselves do not match the work done by Crater, but there is a nearly negligible difference 
between them, which is an expected result and thus some theories assume this ab initio{^). It is also 
interesting to note that the up and down quark masses are almost right in between the masses given by the 
two major sources of inspiration and material for this work. Crater and Isgur f |39[ [5]). 

On the whole though, the fit is nearly as accurate as others, most notably the work of Capstick and Isgur 
([6]), which is generally regarded as one of the more valuable references for theoretical baryon spectroscopy. 
The only marked difference of the results of our model versus other models is that while our higher order 
baryons are generally fit fairly well, the lower ones are not fit as well, relatively speaking. However, as is 
discussed below, this may actually reinforce that the fundamental approach is sound. 
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6.4 Conclusion and Future Work 

The model has shown that it is possible to use the purely two-body approaeh for a good fit of the baryon 
spectrum. However, some work remains to be done though to truly count this model as accurate for the 
purposes of making quark model predictions. This should not bee seen as a detriment to our model but 
rather a triumph of it as we are aware of interactions, such as how the interactions of two particles affects 
the third, that we did not consider and may be significant. The high masses of the octet can largely be 
attributed to our potential model not discriminating between such states as the N and the A with anything 
more than the spin-spin interactions, which may simply not be enough to account for approximately a 
300 MeV difference between baryons. All of the other interactions require either a radial or an orbital 
excitation. 

As for future work, it may be possible to introduce fully three-body forces in addition to the two-body 
ones and to use a fully three-body approach for a coordinate system and JLS couplings. Total JLS 
couplings for a three-body system are usually done in a mathematically rigorous fashion by coupling two 
particles together and then coupling their Clebsch-Gordan coupled two-body system to a third particle for 
a complete three-body system. A fully three-body approach to angular momentum couplings may at the 
very least yield a more elegant formalism and perhaps better overall results. A system derived purely for a 
three-body problem and including three-body JLS couplings may include additional interactions not seen in 
a two-body model. We believe this will solve the issue of the same family of particles (i.e. E, A, N) lacking 
in enough differentiation. Realistically speaking, it can be viewed as positive in this work that the lower 
order baryons are not fit as well as other models, but rather the model as a whole appears to be shifted. 
This lends credence to the theory as a whole being fundamentally sound, but merely incomplete. In other 
words, if the "missing piece" were added to the model in such a way as to affect all the baryons, it woiild 
be more likely to shift all of them in the correct direction, as opposed to if some of them were nearly exact 
in the current fit and thus an overall change/addition to the model would be likely to shift them away from 
values in good agreement with experiment. 

In addition, it may be possible to extrapolate out to an infinitely large basis and thus get an exact value 
for the ground state, rather than simply truncating our basis size at an arbitrarily large value. This may 
yield a better overall fit, but it is worth noting that the truncation was done at the N=3 value because higher 
values did not appear to change substantially, so it may not give any more accurate results. 
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7 Appendix 



A Explicit forms of the Potentials 

The full two-body potential is written as 

where the individual i — j interacting potentials are given by (r being the distance between each pair in that 
particular i — j combination) 



2(J-' + l/r)(cosh2/C-l) ^,2 2/C'sinh2/C 

r r 

r (J'' + l/r)(cosh2/C- 1) /C'sinh2/C 

•Pso — 1 , 

r r r 

q>soD = (Z'cosh2/C-g'sinh2/C), 

*sox = (g' cosh 2/C - /' sinh 2/C), 

2/C'sinh2/C 2(J-' + l/r)(cosh2/C - 1) 2J''/C' V^/C 



3r 3r 3 3 



. + (3^'-^^ + 3/r)sinh2/C ^ (-^' - + lA)(cosh2/C - 1) ^ 2^,^, ^ ^2^j^ 

+ sinh 2/C _ ^^^^ 
r r r 

These in turn are made up of the functionfUl 

n{r) = \v'g + 'J^'^ + {G' - 21C% 

m(r) = -iv2g+^g'2 + g'^'-/C'2, (A3) 



^^These are all two-body potentials acting in a three-body system, so the subscripts 1 or 2 refer to the 1st or 2nd particle in 
that pair, not the first or second particle in the overall system. 
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1 ^2M^^-^lMi 



2 



2 



2 e2'7ii+eim2 

{L' ^g'){E2M2 + EiMi) 

2{E2Mi+EiM2) 
{V^L-V-'g){E2M^+E^M^) _ _ 2 K - mlf 

2{E2Mr+EiM2) 2 (^sMi + iJiMj)^ 

£■2 Ml + E1M2 

E'^Mi + E[M2 + E2M[ + EiM'2 



Mi 
Mi 

Ei 



/ 2 2m^5 + 52 
Y ™' ^ {w-2A)/w 

sjm\ + 2m^5 + 52 

(m^ + 2m^5 + ^2)3/2 Mf 
- A 

- 2A)lw 



^el - 2e^A + 

i-e^A" + A") i-e^A' + A')^ ^ js^A" + A") E 

^ef-2s^A + A^ (£2-2£^A + A2)3/2 £ 

-^Og( (l-2AM ) 
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s' 

A' 



w -2A 
47rA2 



27 

6162 Stt IGBtt 



47rr2 27rHog{e^K + 21rHe^K + ^W{log{e^K + ^)) 



2 



_ 16B7t{-4B + {B + 3e^Kr^A^)log{e^K + ^)) 



171x1712 {w^ — m\ — m^) 

W ZW 

^ _ ^J7n{ + exp(2g)(2m^5'+5'2) + ^Tnj + cxp(2g)(2TO^5+5'2) 

mi + TO2 



MiAf2 V ^i' - 2yl (w - 2 A) 



M1M2 M1M2 
V^5(to^ + 5') + 45"(m^ + + (2m^5' + S'^)V'^ A 4(2to^5 + 5'^)^^^ 

W-2A ^ [w - 2Af ^ (u; - 2Af 



(AlO) 



The scalar and vector invariants A(r) and S(r) form the basis for all of the above functions. They are given 

by 

27r\ogiKe^ + j£^) Anr 



B Spin-Flavor-Space States 

This appendix contains the derivations of the spin- flavor states for all the baryons in our fit. They are 
composed of products of spin wavef unctions, denoted x, and flavor wavefunctions, denoted as </> |47) . The 
spin wavefunctions are, explicitly 







=ttt, 












;n), 
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- m - in)- 


(Bl) 



Table 8: Baryon flavor wavefunctions 







7 




p 




^ — duu) 


-^{2uud — duu — udu) 


n 




{udd — dud) 


{dud — udd — 2ddu) 


A 




{usd + sdu — sud — dsu — 2dus + 2uds) 


^{sud + usd — sdu — dsu) 


A++ 


uuu 






A 1 

A+ 


uud 






aO 


udd 






A — 

A 


J J J 
add 








uus 


f dill 11 Cl l\ 


( CI 1 '} 1 1 1 CI 1 ^ 1 11 1 c\ 






— { QiirJ -i- criii 7/0/7 rlQli\ 


^ ( 11 cfi -\- Qflii — t- Qii rl —1— ri Qti Q ^oirl q\ 




ilUjO 




_L ( aflrl — fla(l — rl d t\ 




USS 


^ ( Q1IQ 


_i_f'Oee'?/ QHQ ')/oo^ 


, — _ 


r/.s-.s- 




_L('9s-s7/ — sy/s — (his) 
Vi'> • ' 












uuc 


-7= fcuu — ucu) 


—7= (cuu — ucu — 2uuc) 

v6 ^ ^ 




uub 


-^{buu — ubu) 


-^{buu — ubu — 2uub) 


Ac 




■^^{ucd + cdu — cud — dcu — 2duc + 2udc) 


^ {cud + ucd — cdu — dcu) 


Afc 




(ubd + bdu — bud — dbu — 2dub + 2udb) 


^ {bud + ubd — bdu — dbu) 




use 


-^{suc — use) 


-^{2scu — sue — use) 


Sb 


usb 


--y=[sub — usb) 


-^{2sbu — sub — usb) 


fib 


ssb 


-^{ssc — scs) 


-^{css — scs — 2ssc) 



Ta.l)l(> 9: Total s})iu-lia\"or-spacc wax cruuctious 



N 


J 


L 


s 


Total State 




8 


i 

2 





i 

2 


^{^'x' + rxl^o 




10 


3 
2 





3 
2 




*2 


8 


i 3 


1 


i 


m<p'x"+ci^"x')i'' + {^'x'-(i>"x")r] 


*3(J=i),*4(J=f) 


8 


2 ' 2 ' 2 


1 


2 




*5(J=i),*6(J=i),*7(^=|) 


10 


1 3 

2 ' 2 


1 


i 

2 


^[(/.^xV' + '/'^x'V"] 


vl/8(J=i),Vl'9(J=§) 


1 


1 3 

2 ' 2 


1 


i 

2 


^rx'V'-rxV"] 


*io(J=^),*n(J=i) 



There are four diflFerent flavor wavefunctions, denoted as <j)" , 
and the singlet state 

t^" = —^{uds + dsM + sud — dus — usd — sdu) (B2) 



There are eleven possible combinations of these spin and flavor states for (most of) the known baryons. 

These wavefunctions then define a grouping of baryons and the individual baryon itself is defined by the 
flavor state from there, as given below 
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Table 10: Baryons and their corresponding spin-flavor wavefunctions 



*i P, N, A, E+, S°, iV(1440), A(1600), S(1660), S(1690), S+(2455), E+, E+, A+, A+(2595), Ag 

*2 ^ A++, A+, A°, A", I]+(1385), I]"(1388), S'(1390), S°(1530), S-(1535), 17", A(1600), 2(1690) 

*3 ^ iV(1535),A(1670),S(1750),E(1880) 

^-4 ^ iV(1520),A(1690),S(1670),S(1820) 

*5 ^ iV(1650),A(1800),S(1750) 

^-e ^ iV(1700),E(1940) 

*7 ^ iV(1675),A(1830),I](1775),S(1950) 

*8 ^ A(1620) 

*9 ^ A(1700) 

^-10 ^ A(1405) 

^-11 ^ A(1520) 



C Operator and 9j Details 



The eleven possible wavefunctions [17] are 



N 


J 


L 


S 


Total State 




8 


1 

2 





1 

2 






10 


3 
2 





3 
2 




*2 


8 


1 3 

2 ' 2 


1 


1 
2 




*3(J=^),*4(J=i) 


8 


13 5 
2 ' 2 ' 2 


1 


3 
2 




*5(J=i),«'6(J=|),*7(J=|) 


10 


1 3 

2 ' 2 


1 


1 

2 




*8(J-^),«'9(J=i) 


1 


1 3 

2 ' 2 


1 


1 

2 




*io(J = 5).*ii(^ = i) 



where "00 ^ ground state Gaussian and ij/ and V'" are orbitally excited Gaussians given by [47^ and each 
have different coefficients for the 9j recoupling (radial excitations are accounted for in the explicit forms of 
the Gaussians themselves). The tables below give the 9j coefficients for every possible coupling that we use. 
They are constructed in a similar manner to a clebsch-gordan table, where we have ( [54] ) 





possible spin for Sjk — > 





1 


possible / for Ij I 













9j coefiicient( Ji, Jjfc) 


9j coefficient ( Ji, Jjfe) 


1 




9j coefficient( Ji, Jjfc) 


9j coefficient ( Ji, Jjfe) 
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where Jj is the total J for the "spectator" particle (i.e. particle one in the 2-3 interaction) and Jjk is the 
total J for the interaction in question, which is what we are really concerned with for working out the state 
explicitly. 



2 2 ^ 





Sjk 





1 















i(io) 




1 




-V2(i,i) + y|(|,i) 


i(i,0) + f (i,l) + i(i,l) + f (1,2) 



J — -, I, — 1, 5 — -, Zj — 1 — Zj 










1 


h 













i(io) 


-i(i,l)-f^/2(|,l) 


1 









1 3 
J = —, L = 1, S = = 1 — Ij 










1 


h 
















-|^/2(i,l) + i(i,l) 


1 







^/i(iO) + ^(i,l) 



J — -, L — 0, — -, — 





Sjk 





1 


h 


















1 







1(1, 0) - iy(2)(i, 1) + ^(|, 1) - iv^(i, 2) + ^(|, 2) 



3 1 
J = -,L = 1,S = = 1 — Ij 
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L) 


1 















1(1,0) 


!(ii) + f (i,i) 


1 




i(ii) 


~75^2' 1) + 3 \/?(2 ' 2) 



3 3 
J = -, L = 1, S = = 1 — Ij 










1 


I, 
















3^2' / 3^2' 


1 







y|(i,l) + ^(i,2) 



5 3 
J = L = 1, S = -,li = 1 — Ij 










1 


I, 
















1(1,1) 


1 







1(^2) 



Since the states are all independent of total M, we ean set M = J for the higher J states in eaeh set of 
wavehmctions. This allows the use of the nmeh simpler ladder operators as opposed to the 9j method, which 
has the side benefit of being a check on both approachcis (many were done both ways to ensure accuracy). 
The states on which this is possible are ^i, ^2, ^'3, *7, *9, and ^'n. The reason it is simple is because any 
state with M = J = \L + S\ will only have the maximum possible value for Mg and so any operator that 
changes total Ms will have a different set of quantum numbers from the state to which it is coupling and 
will thus give a matrix element of zero. The simplified versions of the ladder operators are (where we have 
used said property of M and Ms) 

Oj-Oj — r 'Jiz'^jz: 

SO = L- S = L,_{S,, + Sj^), 
T = sin' + 3Si,Sj,cos^0- Si -Sj, 

SOD = L^Siz-Sjz), 

SOX = L,^ Si-Sj+-Si+Sj. ^_ ^^^^ 
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Note both that the spin-spin operator is simple regardless of the state and thus we can use it either way on 
any state and that the spin-orbit tensor operator does not appear. This is due to the fact that the spin-orbit 
tensor operator does not have a simple form, even with the simplification that any change in total Mg will 
be zero. Using both of these approaches, we can reduce all of our operators to simple numerical factors onto 
the radial wavefunction, which vary depending on the state. The tables are written in the format 



total wavefunction 


interaction 


quark coupling 


numerical coefficient 



So that the table for (for example) would read that all interactions but spin-spin are zero and all spin-spin 
interactions give back a coefficient of 



*1 


SS 


SO 


T 


SOT 


SOD 


SOX 


1'2 


SS 


SO 


T 


SOT 


SOD 


SOX 


12 


1 

4 

















12 


1 

4 

















13 


1 
4 

















13 


1 
4 

















23 


1 
4 

















23 


1 
4 

















^3 


SS 


so 


T 


SOT 


SOD 


SOX 




SS 


so 


T 


SOT 


SOD 


SOX 


12 


1 
4 


4 
3 














12 


1 

4 


2 
3 














13 


1 
4 


4 
3 














13 


1 

4 


2 

3 














23 


1 
4 


2 
3 














23 


1 

4 


2 
3 








1 

2 






^-5 


SS 


SO 


T 


SOT 


SOD 


SOX 




SS 


SO 


T 


SOT 


SOD 


SOX 


12 


1 

4 


10 
3 


-2 


1 








12 


1 

4 


4 
3 


8 
5 


8 
5 








13 


1 

4 


10 
3 


-2 


1 








13 


1 

4 


4 
3 


8 
5 


8 
5 








23 


1 
i 


5 
:S 


-1 










23 


1 
1 


_ 2 


4 


4 










SS 


SO 


T 


SOT 


SOD 


SOX 




SS 


SO 


T 


SOT 


SOD 


SOX 


12 


1 

4 


2 


2 
5 


2 
5 








12 


1 

4 


4 
3 














13 


1 

4 


2 


2 

5 


2 

5 








13 


1 
4 


4 
3 














23 


1 
4 


1 


1 
5 


1 
5 








23 


1 
4 


-1 








1 
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Wg 


bb 


bU 


i 


bUi 


bUU 


bUA 




bb 


bU 


1 


bU i 


bUU 


bUA 


12 


1 

4 


2 
3 














12 


1 

4 


4 
3 














13 


1 

4 


2 
3 














13 


1 

4 


4 
3 














23 


1 

4 


1 

2 








1 

2 





23 


1 

"4 


4 
"3 








-1 






fll 


SS 


SO 


T 


SOT 


SOD 


SOX 


12 


1 

4 


2 
3 














13 


1 
4 


2 
3 














23 


1 
4 


1 
6 








1 

2 






D Coordinate System Detailed Derivations 

This appendix contains explicit details for the coordinate transform derivation. First, we define 
relative coordinates 

p = r2 - r3 

We2 , ^£3 

^ = 7 \ \~''2 + 7 ; r— ra 

(£2 + S3)Sl (£2 + S3)Sl 

ri - Fa = ■ p - A 

£2 +£3 

ri-ra = ■ p + \ 

£2 + £3 

r2 - rs = p 

where the e's are the energies of each quark and 

W = Si + S2 + £3, 



as usual. Our base wavefunction is 



= ]^plp)lxg^npalp^/2-nxal\^/2Y^PY^>^ 



97 



which must be normahzed. As usual, the spherical harmonics are already normalized 



' Y 47r (/ + m)! ' ^ ' 

SO we just have to normalize the radial part of the wavefunction, which gives 

,2 r[(2/p + 3)/2] r[(2/A + 3)/2] 



= N 

(npa2)(2;p+3)/2 („^^2)(2i,+3)/2 

[hi, 

N = 



(npa2)(2ip+3)/2(^^Q,2)(2i^+3)/2 



T[{2lp + 3)/2]T[{2lx + 3)/2] 



I (n.a2) (2ip+3)/2 ( 2 )(2/^+3)/2 ^ 

" V r[(2Zp + 3)/2]r[(2ZA + 3)/2] ^ 
or more completely with the appropriate Clebsch-Gordan coefficient 



r[(2Zp + 3)/2]r[(2/A + 3)/2] 



mpmx 
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This wavefunction is now sufBcient to perform the r23 = p integration and the kinetic terms (which are 
analytic). The kinetic terms are given by 



dp 



V2] 



(*p|T|*p') 



(*,|T|V) = 



l,-2^-in,+n')c.lp^/2 



p'-+^N,N2[lp{lp + l)p'^ 



-{2lp + 3)npa>'-e-("-+"p)">'/2 ^ „2^4^'p+2e-(".+"P">V2] 



/(npa2)(2/p+3)/2„^(2Zx+3)/2 / („/^^2 ) (2;p+3)/2„'^(2'A +3)/2 

/ r[(2;p + 3)/2] V r[(2/p + 3)/2] 



[lp{lp + 1) 



r[(2Zp + 3)/2] ^ , , r[(2;p + 5)/2] 



T[{2lp + l)/2] 



'2(K + nyQ:2/2)(2'p+3)/2 P ^2([np + n;,]a2/2)(2ip+5)/2 
X J2{lphmpmx\LM)Y^'^^ x ^(/;;»'jLM)r7^ 



2([r7.p + n;,]a2/2)(2ip+i)/2 
(D8) 



(D9) 



In order to perform the ri2 and ri3 integrations, we redefine the wavefunction in what we call spherical 
tensor form as 

= iVp'pA'^e~"''"p''^/^""^"^'^^/^ 

X {lplxmpmx\LM)Y,^%^- 

= ^e-^.o.yi2-n,o.lx-l2 ^ ^l^l^mpmx\LM)p^\^^ (DIO) 



where 



(Dll) 



For ri2 integration, we will define 



ri2 = ri-r2=p' (D12) 

, wei we2 
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in order to create an integration variable purely in terms of ri2 = p' ■ Then we rewrite our wavefunction in 
terms of new variables and as tensors 



eiri +£:2r2 £2 + ^3 , £1 „^ 
r2 - rs = r2 = vi ^ ri (D13) 

£3 £3 £3 

(£2 + £3)£i £i£2 N / _|_ , £2 + £3 £1);^' 

£3(£i + £2) £3(£i+£2) £ £ 



£1 / , \' 
P = ^ P 

£1 + £2 



^ _ W£2 ^ ^ W£3 ^ _ W£2 ^ ^ W£3 -£iri - £2r2 , (D14) 

£l(£2+£3) ^ £l(£2+£3) ^ £l(£2+£3) ^ £l(£2+£3) £3 



W , W£2 W£2 s 

-ri + (— ■ ■ -)r2 



£2 +£3 £i(£2 + £3) £i(£2+£3)' 

, W£2 , £3 

A = — -p A 

(£2 + £3)(£l + £2) £2 + £3 

A = W£2 , £3 y, 

(£2+£3)(£l+£2)^"^ £2+£3 

A general expectation value of any potential $12 which depends only on the coordinate ri2 is 

(*„|$12(ri2)|*™> - (*„|$12(p')l*m> (D15) 

= iVi7V2 / (^pi™p"^AiiAf)p:„^A;,^ a;i™XAi^'^^>™,A™,e-^("^+";')">'-^("^+"'^)"'^^'a>i2(p')dVrf'A 

(D16) 

Using above identities for p and A (jD13ID14ID12p we arrive at 

(^'„|*12(p')l*m) 



X 1: (,-,,„■;„.- |L-M-)(^A... + A-„.j( A.,. - ^a:,,j (dit) 

m'm' 
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where for simplicity 



2 'ei+£2' 2 '(£2+£3)(£i+£2) 

, v-p + "'p)"L £1 . , , /^ 2 w;£2£3 



{Hp + n')al ex 

2 ^7r^.> ^ ^ '"^^"^ (£2 + £3)^(£l + £2) 



^ ^ K + ^ K + n^Jaj ^ £3 .2 pi8) 

2 2 £2 + £3 

In order to eliminate the cross term in the Gaussian, we make one final change of variables and essentially 
complete the square, so a variable x is defined as 

X = A' + Ap' (D19) 

w _ _ 1 , 

Prom this point on the calculation will continue for Z = 1 as it becomes extremely complicated for general Z, 
but the formalism is the same. The expectation value now is 

(*„|$12(P')I*™) 



m' m'. 



mpmx 

. W£2 I £z , b I £1 , ^ / \* 

''^(£2+£3)(£l+£2)^'"^ ~ i^Ti^^"''"^ ~ 2^^™^''^ ^^7^1^^-^+^™" - Yc^rnJ 

WS2 , ^3 / ^ 



m'm\ 
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Expanding this out, we obtain 



($12) = N1N2 J e-<^'+^^p'^"+^''-'&^p'"d''p'd^x' {UmpmxlLM) 



nipmx 

^(£2 + £3)(£i + £2)^"^ £1 + £2^"" (£2 + £3)(£i + £2)^™-''"^ (£2 + £3)(£i + £2)^"^ 2c' 

^3 * ^1 ^3 if ^ , ^3 



£2 + £3 "^"^ £1 + £2 £2 + £3 "^"^ ^ £2 + £3 "^"^ 2c ^" 



-I- I * "I /"TiOO^ 

"2c '' El + £2 2c ~ ^^^'"A/'nV-l 



X 



W£2 , £1 , , W£2 , We2 

-P ' P • + 7 77 7P ' 



'(£2 + £3)(£l+£2) "A£l+£2 (£2 + £3)(£l+£2) "a (£3 + £3)(£i + £2) "a 2c^"V 

£3 £1 / £3 , ^3 b 



£2 + £3 "'^£i+£2""- £2 +£3 ^ ^ £2 +£3 "'^2c^™p 

Now that we have expanded it we wiU go backwards and replace the spherical tensor form with the explicit 
forms of spherical harmonics from 13021 and the coordinates 



($12) = iViiV2ye-^(^'+*''')'+("-^)''"dVrfV ^ (llmpmAiLM) 



r™^*(p')n"'*(p>" + 7 T7^ -Yr'-*{p)Yr'*{x)p'x 



(£2 + £3)(£l + £2) £1 + £2 (£2 + £3)(£l + £2) 



(£2 + £3)(£l + £2) 2c £2 + £3£l+£2" 



£2 + £3 £2 + £3 2c 

+ + )i^r*(p )P'^](D24) 
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X (1Z»'JL'M') 



'■p^"■\ 



(£2 +e3)(ei +£2) ei +£2 '■'^ ^ ' ^ (£2 + £3)(£i + £2) " 

if£2 ■xa"»1*/"/\v™p* /"/^ 12 ^3 £1 



(£2 + £3)(£l +£2) 2c £2 + £3£l+£2 

^3 y-l*(^)y->(^)a.2 + ^^^y-l*(^)y->(^')^'a. 



£2 + £3 £2 + £3 2c 

+ lYr'^*cp')Yr'^*ix)p'x ^^Yr'^*cp')Yr''*cp)p"] m 



We use the following identities of spherical harmonics to perform the angular integrations and simplify 



Yp{e,cf>)YJ^%e,<P) = E V ^^"^4^(2^+1)^ {jihm,m2\jij2jm){j,j200\jO)Yp{e, <f>) (D26) 



J (e, <f>)Y,^^ {6, cP)Y,r (0, <t>)d9dcl> = sf- 





h 


:( 










roi m2 m3 

(D27) 

So now 

($12) = A^iA^2y" e-'=(^)'+("-^)^"d^a;dV E (11™,.™^!^'^) ^ (llm;m';^|I,'M') (D28) 



TTlpTTlx 1 1 1. ^ 1 1 1, 



-n^X*/;,/^^m^/V^V'"''*/'?,'^^'"p/5./^„/4 (-2c£l + &(£2 + £3))^(K£l + £2)£3 + 2ce2W)^ 



16c4(ei+£2)2(£2+£3)^ 

+p-x^Yr*{i)Y< {£)Yr^*{p')Y< ip') 

_l_^2 ,2 £3(-2c£i + 6(£2 +£3))(^(£l + £2)£3 + 2c£2w) 
c2(£l +£2)(£2 + £3)^ 



X [y/"^* (f )yi"''* (f )y/"^ (pO^i"" (pO + ^r'* (pO^r" (p')^r' (^)n'"''* (^) 

+y/"^* (x)yr^ {x)Yf^ (p')yr'* (p ) + yr^* (p )i"r''* (p )^r^ (^)^r'' (^)] 

+x2p'2y™- (p')yr {p')Yr*{x)Yr'' (i) ^ ^ 



4c2(£i+£2)2(£2+£3) 



2 



2 



+x^Yr\x)Yr'' imr^* imr- ) (^^f^ ) 

This is a gigantic equation when everything is worked out, so we will treat it term by term. The first term 
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gives 



(llm,m,|LM) J] (llmXAl^'^'^On"^*!/'' ) 

,4 (-2c£i + b{e2 + £3))^(&(£l + £2)g3 + 2c£2w f 

16c4(£i+e2)2(e2+£3r 
^ (llm,m,|iAf) 5] (ll™X,|L'Af')n™^*(p)n™^*(p) 



mpfiix 



(-2c£i + 6(£2 + £3))'(&(£i + £2)63 + 2ce2wY 

16c4(£i+ £2)2(62 +£3)4 ^ ^ 



Again from Eq |D281 we have 



^ (llmpmAiLM) ^ (llm;m';,|L'M')x4y™^*(x)yi'"^(x)yr''*(f)rr''(f) — 
/ (£2 



mpmx 



£3) 



2 



2 



4 ^3 



(£2 +£3)2 



J2 {llmpmxlLM) ^ (llm;TO';^|L'M')y™''*(x)y™^*(a 



5^ V ^(2; + 1) (ll'^W;illM(1100|11^0)>r(P ) (D30) 



So the expectation value is 

($12) = N1N2 J e~''^'''^'+^''~'r^^p'^d^xd^p' J2 (ll"ip™A|iM) J2 {llm'pm\\L' M') 



mp-mx 



^ ,4 (-2e£i + b(£2 + esffibjei + £2)£3 + 2c£2w)2 

16c4(£i+£2)2(£2+£3)4 

+p-.2yr^* (i)^'"'^ (-)n™^* ip)yr'' ^p) '"^~%tTe,V'^^' 

, t 2 12 , 2 12 , 2 12 , 2 I2\^'^ ""^^ 
+ {X p + X p + X p + X p ) 



2 ,2 , 2 /2 , 2 /2 , 2 ,2^£3(-2c£i+6(£2+£3))(fo(£l+£2)£3 + 2c£2w) 



4c2(£i +£2)(£2 +£3)3 
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X [y™^* (f )y"''* {x)Y{"^ {p')Y^' ip') + y;"^*{p')y^' (p')^r' {x)y^''* {x) 

+Y{^'*{x)Y^'' ip')Y;"'*{p') + (P )^r^ {x)Y^'' {x}] 

wp''Yr*{p')Yr'^ {p')Yr*{x)Yr'" (*) ^ ^ 



Im 



We can now perform the angular integrations for the spherical harmonics, which are identical for both 
the p'* and x'^ term as there is no dependence on the angular pieces in the potentials 



mpmx m'^m'^ Im • 

x(llm';^m'|llZm)(1100|llZ0)y;'"(p')) 



{21 + 1) 



mpmx m'^m'^ Im y 

( 



X 



2(1 + 1)(21 + 1)(2Z + 1) 
47r 



11/1 1 Z 
/ 1 m\ nip m 



(D32) 



All other terms either integrate to zero (due to an odd number of Yjj„'s) or to Kronecker deltas, so the rest 
of the integral is just a polynomial, given by 

($12) = NiN2 j e-'=(^)'+("-^)''"a;V2d:cdp' ^ {llmpmx\LM) ^ {llm'pm'^\L' M') 



mpmx 



X 

Im 



X 



2(1 + 1)(21 + 1)(2/ + 1) 1 1 M 1 1 / 
\ j \ mx mp m 

(-2c£i + b{e2 + £3))2(6(£i + £2)53 + 2C£2W)2 
16c4(£i +£2)2(£2+£3)4 

,2 2ei(-2c£l + 6(£2 +£3))2 



4c2(£2 + £3)4 



_^^^2 /2 £3(-2c£i + 6(£2 + £3))(&(£l + £2)53 + 2c£2W) 
4c2(£i +£2)(£2 +£3)3 
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, 2 12. r (&(gl + g2)g3 + 2cg2w)^ 

A A P P 4c^(£-^ +£2)^(62 +£3) 



2(1 + 1)(21 + 1)(2; + 1) 
47r 



1 1 I 




1 1 I 

m\ rup m 



(D33) 



Performing the x integration yields 



($i2> = ^1^2 / *i2e-'=(^^'+("-^)'''V"rfp' {llmpmx\LM) ^ {\lm'pm'^\L' M' 



){ 



2(1 + 1)(2H-1)(2/ + 1) 



4-77 



1 1 ? 





1 1 I 



m\ nip m 



xp 



,4 (-2c£i + 5(£2 + £3))2(6(£i + £2)£3 + 2c£2w)^ ^ 



+P'^X 



16c4(£i +£2)2(£2+£3)^ 

,2 2ei(-2c£i + 6(£2+£3))'3 ^ 1 



4c3/2 



4c2(£2+£3)^ 4"^ 2c5/2 



/2 /:r 1 £3(-2c£i + 6(£2 + £3))(6(£i +£2)£3 + 2C£2W) 



2c5/2 



4c2(£i+£2)(£2+£3)3 



"I" 2^to^^j^(5to^TOp "I" ^m'^m'p^mxmp] 

(&(£l + £2)£3 + 2C£2W)^ 



4c2(ei+£2)2(£2+£3)2 



1+1+1+1 



(£2 + £3)' j 



/ 2(i + i)(2r, + i) 

47r(2/ + 1) 



{ll'^m'^m' \ll'^lm){ll')^00\niO) 



2(1 + 1)(21 + 1)(2/ + 1) 
4tt 




(D34) 



So now the expectation value is in terms of just one variable, which can be numerically integrated easily. 
Doing the 13 just involves changing the a,b,c constants and the epsilons and the 23 is trivial since it 
doesn't require any additional substitutions. 

We can follow a similar method for Z = 2, by constructing new states as functions of the I = 1 states. 
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The explicit forms of the spherical harmonics for ^ = and / = 1 are 




(D35) 



which can be used to rewrite the I — 2 spherical harmonics as 



y-2 



^2 



^^2 



sin6lcos6le^''^ = Y°Y^-^2<^, 



1 15 
iV 2^ 

1 /T5 

2 V 2^ 



{Y,yv5^, 



1 15 
~2V 2^ 
1 [15 
iV 2^ 



9cos0e*^ = ri"yii2y^, 



So then, since our wavefunction for Ip — Ix — 2 is of the form 



l^n(P,A)) = 



15i/ TT 



X {22mp7n,x\LM)Y^''Y^\ 

niprnx 

we can use the original tensor substitution of Eq. (|Dlip to define 



P'Y,-' 
p'Y,-' 



9 ^2 



P ^2 



P "2 



3 



V3' 

P^rW2/|.p„p,2/|, 



V3VT' 



(D36) 



(D37) 
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Therefore, we can follow the same procedure outlined for I = 1, the only major difference being that the 
algebra is more complicated. However, since we're going to have to do at least one numerical integration 
anyway, we can use a program such as Mathematica (or Wolframalpha) to analytically integrate what is 
analytically integrable and then numerically integrate the rest without doing an absurd amount of algebra 
(this is what was done in this work). This method should actually hold to allow for higher order I states to 
be done relatively simply, since essentially all that will become more complicated is the algebra arising from 
the substitutions (this may be part of future work on this subject). 
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